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Abstract 



In the framework of the rest-frame instant form of tetrad gravity, where 
the Hamiltonian is the weak ADM energy Eadm-, we define a special com- 
pletely fixed 3-orthogonal Hamiltonian gauge, corresponding to a choice of 
non-harmonic 4-coordinates, in which the independent degrees of freedom 
of the gravitational field are described by two pairs of canonically conjugate 
Dirac observables (DO) ra(r, <?), 7ra(r, (t), a = 1,2. We define a Hamiltonian 
linearization of the theory, i.e. gravitational waves, without introducing any 
background 4-fnetric, by retaining only the linear terms in the DO's in the 
super-hamiltonian constraint (the Lichnerowicz equation for the conformal 
factor of the 3-metric) and the quadratic terms in the DO's in Eadm- We 
solve all the constraints of the linearized theory: this amounts to work in a 
well defined post-Minkowskian Christodoulou-Klainermann space-time. The 
Hamilton equations imply the wave equation for the DO's ra(r, a), which re- 
place the two polarizations of the TT harmonic gauge, and that linearized 
Einstein's equations are satisfied . Finally we study the geodesic equation, 
both for time-like and null geodesies, and the geodesic deviation equation. 
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Typeset using REVT^ 



I. INTRODUCTION 



In a series of papers [1-3] the rest-frame instant form of the Hamiltonian ADM formu- 
lation of both metric and tetrad gravity was given. The aim of these papers was to arrive 
at a formulation of general relativity with matter such that the switching off of the New- 
ton constant G would produce the description of the same matter in the rest-frame instant 
form of parametrized Minkowski theories [4,5,1], with the general relativistic covariance de- 
parametrizing to the special relativistic one of these theories, and to study its Hamiltonian 
formulation consistently with Dirac-Bergmann theory of constraints. 

To implement this program we must select the following family of non-compact space- 
times: 

i) globally hyperbolic, so that the ADM Hamiltonian formulation [6] is well defined if we 
start from the ADM action; 

ii) topologically trivial, so that they can be foliated with space-like hyper-surfaces 
diffeomorphic to (3+1 splitting of space-time with r, the scalar parameter labeling the 
leaves, mathematical time); 

iii) asymptotically flat at spatial infinity and with boundary conditions at spatial infinity 
independent from the direction, so that the spi group of asymptotic symmetries is reduced 
to the Poincare' group with the ADM Poincare' charges as generators ^. In this way we 
can eliminate the super-translations, namely the obstruction to define angular momentum 
in general relativity, and we have the same type of boundary conditions which are needed 
to get well defined non-Abelian charges in Yang-Mills theory, opening the possibility of a 
unified description of the four interactions with all the fields belonging to same function 
space [5] . All these requirements imply that the allowed foliations of space-time must have 
the space-like hyper-planes tending in a direction-independent way to Minkowski space- 
like hyper-planes at spatial infinity, which moreover must be orthogonal there to the ADM 
4-momentum. But these are the conditions satisfied by the singularity-free Christodoulou- 
Klainermann space-times [7], in which, in presence of matter, the allowed hyper-surfaces 
define the rest frame of the universe and naturally become the Wigner hyper-planes of the 
rest-frame instant form of the parametrized Minkowski theories describing the same matter 
when G ^ 0. Therefore there are preferred asymptotic inertial observers, which can be 
identified with the fixed stars. These allowed hyper-surfaces are called Wigner- S en- Witten 
(WSW) hyper-surfaces, because it can be shown that the Prauendiener re- formulation [8] of 
Sen- Witten equations [9] for triads allows (after the restriction to the solutions of Einstein's 
equations) to transport the asymptotic tetrads of the asymptotic inertial observers in each 
point of the hyper-surface, generating a local compass of inertia to be used to define rotations 
with respect to the fixed stars ^. Besides the existence of a realization of the Poincare' group. 



^When we switch off the Newton constant G, the ADM Poincare' charges, expressed in S^-adapted 
4-coordinates, become the generators of the internal Poincare' group of parametrized Minkowski 
theories [1]. 

^Instead the standard Fermi- Walker transport of the tetrads of a time-like observer is a standard 
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the only other needed property is that the admissible WSW hyper-surfaces must admit an 
involution [10] allowing the definition of a generalized Fourier transform with its associated 
concepts of positive and negative energy, so to avoid the claimed impossibility to define 
particles in curved space-times [11]. 

iv) All the fields have to belong to suitable weighted Sobolev spaces so that the allowed 
space-like hyper-surfaces are Riemannian 3-manifolds without Killing vectors: in this way we 
avoid the analogue of the Gribov ambiguity in general relativity and we can get a unification 
of the function spaces of gravity and particle physics. 

After all these preliminaries it is possible to study the Hamiltonian formulation of both 
ADM metric [1] and tetrad [2,3] gravity ^' with their (8 and 14 respectively) first class 
constraints as generators of the Hamiltonian gauge transformations. Then it is possible to 
look, at least at a heuristic level, for Shanmugadhasan canonical transformations [12,4,5] 
implementing the separation between the gauge variables and the Dirac observables (DO) for 
the gravitational field ^. A complete exposition of these topics is in Refs. [1-3], where it is 
shown that it is possible to define a rest frame instant form of gravity in which the effective 
Hamiltonian for the evolution is the weak ADM energy Eadm ^ [17]. 

Let us consider tetrad gravity. In Refs. [2,3] there is a new parametrization of tetrad 
gravity, still utilizing the ADM action, which allows to solve 13 of the 14 first class con- 
straints. After an allowed 3-1-1 splitting of space-time with space-like hyper-surfaces S^, we 
introduce adapted coordinates ^. The arbitrary cotetrads appearing in the 4-metric 



of non-rotation with respect to a local observer in free fall. 

^More natural for the coupling to the fermions of the standard model of particles. Moreover 
tetrad gravity is naturally a theory of time-like accelerated observers, generalizing the ones of 
parametrized Minkowski theories. 

^See Refs. [13-15] for the interpretation of the gauge variables as generalized inertial effects and 
of the DO's (the non-local deterministically predictable physical degrees of freedom of the gravi- 
tational field) as generalized tidal effects. The non-locality of DO's (all 3-space is needed for their 
determination) may be interpreted as a form of Mach principle, even if our boundary conditions are 
not compatible with the standard interpretation of such principle given by Einstein and Wheeler, 
who choose spatiaUy compact space-times without boundary [16]. 

^Therefore the formulations with a frozen reduced phase space are avoided. The super-hamiltonian 
constraint generates normal deformations of the space-Uke hyper-surfaces, which are not inter- 
preted as a time evolution (like in the Wheeler-DeWitt approach) but as the Hamiltonian gauge 
transformations ensuring that the description of gravity is independent from the 3-1-1 splittings of 
space-time like it happens in parametrized Minkowski theories. 

^Instead of local coordinates for M*^, we use local coordinates on i? x S « with 
Ta K. B? [x^ = z'^{a) with inverse = a^{x)], i.e. a Y^r-adapted holonomic coordinate basis 
for vector fields Ba = ^ e T{R x S) ^ b'X{(7)d^ = ^£^9^ G TM"^, and for differential 



4 



of the ADM action principle, are rewritten in adapted coordinates, ^E^\ and replaced at 
the Hamiltonian level by the lapse N and shift N(^a) — ^^{a)r functions, cotriads ^e(^a)r ^-nd 
boost parameters ip^"'^ 
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one-forms dxi^ G T*M^ ^ da^ = b^{a)dxf' = ^^^dxf" G T*{R x S). The 4-metric has Lorentz 
signature e (+,—,—, —) with e = ±1 according to particle physics and general relativity conventions 
respectively. The induced 4-metric and inverse 4-metric become in the new basis 



'9AB = Vgrr = e{N^ " ^ grsN' Nj,^ g^r = -e^grsN';^grs = -e^grs} = 
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and Ia = NOat = N5^^ = (7V;0). We 



For the unit normals to S,- we have I 

introduce the 3-metric of S,-: ^g^s = —e^grs with signature (++-I-). If is the inverse of 
the spatial part of the 4-metric (^7™^5us = ^Di the inverse of the 3-metric is ^g^^ = —e^Y^ 



^igTu Sg^^ = SI). We have the following form for the line element in M 



4. 



ds'^ = e{N^ - ^grsN''N'){dTf - 2e^ grsN' drda' - e^grsda'da' 



= e 



N\dTf - ^grsida' + N'dT){da' + N'dr) 
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The conjugate momenta are tt ^ (r, ct) , tt^^^ (t, ct) , ^tt ("^^ (r, ct) , tt^^^^ (t, a) , respectively. There 
are ten primary constraints 7f^(r, ct) ^ 0, 7t^-^{r,a) ^ 0, 7r(^)(T, ct) f=:i 0, ^M(a)(T, ct) = 
^(a){6){c) ^e(fe)r(r, a) ^ttJ^^^{t, ct) ~ and four secondary ones: the super-hamiltonian constraint 
7i(r, (t) ^ and the super-momentum constraints 7i^(r, c?) ^ 0. Therefore there are 14 
arbitrary gauge variables, four of which arc the lapse and shift functions. All the constraints 
are first class and the Dirac Hamiltonian depends on 10 Dirac multipliers. 

As shown in Ref. [1], a consistent treatment of the boundary conditions at spatial infinity 
requires the explicit separation of the asymptotic part of the lapse and shift functions from 
their bulk part: N{T,a) = N(^as){T,^) +'>^{t,^), Nr{T,a) = N(^as)r{T,^) +^r(T, c?), with 
n(r, a) and nr(r, a) tending to zero at spatial infinity in a direction- independent way. On 
the contrary, N^as){T, ^) = -K{t) - \ Ku{r) cr^ and N(^as)r{r, B) = -Xr{T) - \ Xmir) cr". The 
Christodoulou-Klainermann space-times [7] , with their rest-frame condition of zero ADM 3- 
momentum and absence of super-translations, are singled out by these considerations. As 
already said the allowed foliations of these space-times tend asymptotically to Minkowski 
hyper-planes in a direction-independent way and are asymptotically orthogonal to the ADM 
four-momentum. The leaves are the WSW hyper-surfaces. They have A^(as)(T, ct) = — e, 
A^(as)r(r, c?) = 0. As in Refs. [1-3], from now on the lapse and shift functions N and Nj. will 
be replaced by —e-\-n and n^. 

It is pointed out in the papers [1,3] that in order to have well defined asymptotic weak 
and strong ADM Poincare ' charges (generators of the asymptotic Poincare' group) all fields 
must have a suitable direction-independent limit at spatial infinity. Recall that the strong 
ADM energy is the flux through the surface at spatial inflnity of a function of the 3-metric 
only, and it is weakly equal to the weak ADM energy (volume form) which contains all the 
dependence on the ADM momenta. This implies [1] that the super-hamiltonian constraint 
must be interpreted as the equation {Lichnerowicz equation) that uniquely determines the 
conformal factor (j){T,a) = {det^g{T, a)y^^'^ of the 3-metric as a functional of the other 
variables. This means that the associated gauge variable is the canonical momentum 7r<^(r, a) 
conjugate to the conformal factor, it carries information about the extrinsic curvature of E,-. 
It is just this variable, and not York's time, which parametrizes the normal deformation 
of the embeddable space-like hypcr-surfaces S,-. Since different forms of 11^ correspond to 
different 3+1 splittings of M"^, the gauge transformations generated by the super-hamiltonian 
constraint correspond to the transition from an allowed 3-1-1 splitting to another one: this is 
the gauge orbit in the phase space over super-space. The theory is therefore independent of 
the choice of the 3+1 splitting like in parametrized Minkowski theories. As a matter of fact, 
a gauge fixing for the super-hamiltonian constraint is a choice of a particular 3+1 splitting 
and this is done by fixing the momentum 7r0(r, a) conjugate to the conformal factor. This 
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shows the dominant role of the conformal 3-geometries in the determination of the physical 
degrees of freedom, just as in the Lichnerowicz-York conformal approach [18,19,16]. 



As explained in Refs. [1,3], following a suggestion of Dirac [20], we restrict our space- 
times to admit asymptotic Minkowski Cartesian coordinates, namely the admissible 4- 
coordinate systems = z^{a'^) have the property x'^ — > ^'^^^ ^) at spatial infinity 

with -2(^^)(t, a) = 2:^(^)(t)+&(^),.(t) cr"^ [(/x) are flat indices]. Therefore the leaves of the admis- 
sible 3+1 splittings of our space-times tend to Minkowski space-like hyper-planes at spatial 
infinity in a direction-independent way. While a:|^-)(r) denotes the world-line of a point 

(centroid) arbitrarily chosen as origin of the 3-coordinate systems on each S^-, the &(^)^(t)'s 

are fiat asymptotic tetrads with equal to the normal to the asymptotic Minkowski 
hyper-planes. To force the existence of these asymptotic coordinates, we must add to tetrad 
gravity ten extra configuration degrees of freedom, namely the centroid 2;[^)(t) and the six 

independent components of the tetrads &(^)yi(T), and of the conjugate ten momenta, p^^^ and 

a spin tensor S^^^^'^\ As shown in Refs. [1,3], this increase of variables is balanced by adding 

ten extra first class constraints determining the ten extra momenta: = pfoo) ~Padm ~ 0; 

X^-^ = J(^) — J ADM ~ 0> where the weak (volume form) ADM Poincare' charges Padm^ 

JaSm and the quantities pf^^ = ^foo)(M) ^'S)' = Koo){,^)^{oo)iu) Sl^)"^ are expressed in 
E,--adapted 4-coordinates. 

In this way the ten extra configuration degrees of freedom become arbitrary gauge vari- 
ables. After the splitting of the lapse and shift functions in the asymptotic and bulk 
parts, this approach [1,3] implies the replacement of the canonical Hamiltonian / d^a nTi.+ 

UrW {T,a) + Xa{t) P^DM + ^ ^ab{t) Jadm the Hamiltonian / d^a nH + Uj-H^ + 

An TTn + Xfir T^H + )"(a) ^(a) + P{a) VT^)] (t, g) + \a{t) [pf^) - P^dm] + \ Aab(t) [ J(^) - J^Em] ■ 

The following boundary conditions (compatible with Christodoulou-Klainermann space- 
times) ensure that the Hamiltonian gauge transformations preserve these asymptotic prop- 
erties defining our class of space-times 

^e(a)r(r, a) ^r^oo (1 + ^)^{a)r + Oir-^'"^), 

^ '^9rs{r, a) = [^e(a)r ^e(„),](r, a) ^r^oc (1 + y)5r. + 0{r-^/^), 

n{T,a)^,^^ 0(r-(2+^)), e > 0, 
nr{T,a) ^r^oo 0(r"'), 
a^a)ir,a)^r^^ 0(r-(i+^)), 

(^(„)(T,a)^,^oo 0(r-(^+^)). (1.2) 



As shown in Refs. [1,3], to get the rest-frame instant form of tetrad gravity, with its 
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WSW hyper-surfaces as leaves of the admissible 3+1 splittings (when p^^^ is time-hke) the 
gauge fixings must be added with the following procedure : 

i) Add three gauge fixings on the boost parameters, namely <^*^"^(r, a) ~ 0: in this way 
we choose cotetrads adapted to E^. The time constancy of this gauge fixing determines the 
3 Dirac multipliers p(^a) {t-, ^) ■ 

ii) Add three gauge fixings, «(a)(T, (?) ^ 0, to the rotation constraints ^M(a)(r, a) ^ 0: in 
this way we fix a reference orientation of the cotriads. The time constancy of these gauge 
fixings determine the 3 Dirac multipliers /X(a)(T, a). 

iii) Add three gauge fixings Xri'T, ^) ~ to the secondary super- momentum constraints: 
this amounts to a choice of 3-coordinates on E,-. The requirement of time constancy of 
the constraints ^^(r, a) ^ will generate three gauge fixings Lpj.{r,a) ^ for the primary 
constraints 7r^(r, a) ~ 0, which determine the shift functions n,f.{T^a) (and therefore the 
coordinate-dependent gravito-magnetic aspects and the eventual anisotropy of light propa- 
gation). The time constancy of the </?r(T, (?)'s will determine the 3 Dirac multipliers Xnr{T-i 

iv) Add a gauge fixing x(t, a) ~ to the secondary super-hamiltonian constraint, which 
determines the form of the space-like hyper-surface S^. Its time constancy produces the 
gauge fixing ipij, a) ~ for the primary constraint 7^n{T, a) ~ 0, which determines the lapse 
function n(r, a), i.e. how the surfaces S,- are packed in the foliation. Now the 3-1-1 splitting 
of space-time is completely determined and the time constancy of ip{T, a) ~ determines 
the last Dirac multiplier A„(t, a). A posteriori after having solved the Hamilton equations 
one could find the embedding z'^{t, a) of the WSW hyper-surfaces into the space-time. 

v) Add 6 suitable gauge fixings on the gauge variables &(^)^(t), so that, after having 
gone to Dirac brackets, we get J^-^ = J^Em \ab{t) = 0. As a consequence of this gauge 
fixing, the gauge variables a^||^)(r) must be replaced by the canonical non-covariant 4-center 
of mass (r) of the universe. All this implies that the asymptotic Minkowski space-like 
hyper-planes become orthogonal to the weak ADM 4-momcntum Padm ^bat the canon- 
ical Hamiltonian becomes / d^a nH + rirH'^ + Xn'n'n + Xnr'^fi + IJ'(a) ^(a) + P(a) T^'(a) ^) ~ 
At(t) [e(oo) - Padm] + ^r{r) Padm with e(oo) = -e ^/^pf^)- Namely only the four first class 

constraints e(oo) — Padm ~ ^adm ^ With the constraints Padm ~ being 

the rest frame condition for the universe. 

vi) Add the gauge fixing r — T(^oo) ~ 0, implying A^(r) = and identifying the mathe- 
matical time T with the rest frame time ^(oo) = P(oo)(m) ^\^)/Hoo) = P{oo){fi) x[^)/Hoo) of the 
universe. The canonical 4-center of mass x^^^ becomes a decoupled point particle clock. 

vii) Either add the natural gauge fixings Jadm ~ to the constraints Padm ~ or work 
by choosing Ar(T) = as a pre-gauge-fixing (see the discussion at the beginning of Section 
V). 

At this stage the canonical reduction is completed by going to Dirac brackets and it can 
be shown [1,3] that the Dirac Hamiltonian for the rest-frame instant form of tetrad gravity 
is the weak (volume form) ADM energy: = Eadm = ~^Padm- becomes the effective 
Hamiltonian for the gauge invariant observables parametrizing the reduced phase space in 
the rest-frame instant form. 
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To find a canonical basis of Dirac observables for the gravitational field in ab- 
sence of known solutions of the super-hamiltonian constraint, we can perform a quasi- 
Shanmugadhasan canonical transformation adapted to only 13 of the constraints and utilize 
the information that this constraint has to be interpreted as the Lichnerowicz equation for 
the conformal factor 0(t, a) = 7^/^^(t, a) = {det^g{T, a))^!^ ^ g^^^'^^/^ of the 3-metric. The 
result of this point canonical transformation is (a = 1, 2) 
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where a(a)(r, a) are angles and ^^(j^a) are a parametrization of the group manifold of the 
passive 3-diffeomorphisms of S,-, describing its changes of 3-coordinates ^. The entries ~ 
are the new momenta corresponding to 13 Abelianized first class constraints: besides the 
lapse, shift and boost momenta, there are three Abelianized rotation constraints tt (r, a) ~ 

and three Abelianized super- momentum constraints 7f^(T, a) ~ 0. 

The Hamiltonian gauge variables are the 13 configuration variables (^("^(r, a), n{T,a), 
nr{T,a), Q;(a)(T, ct), ^^{T,a) (they depend on the cotetrads and its space gradients) and 
the momentum 7rf^(r, a) conjugate to the conformal factor (it depends also on the time 
derivative of the cotetrads). The variables ^^{t. a) and 7r0(r, a) can be thought as a possible 
4-coordinate system with the Lorentz signature given by the pattern 3 configuration + 1 
momentum variables. 

Even if we do not know the expression of the final variables «(„), tt^, r^, tTq in terms 
of the original variables, the point nature of the canonical transformation allows to write 
the following inverse relations (the form of the cotriad was determined by solving the multi- 
temporal equations for the gauge transformations generated by the first class constraints 
[3]) 

(r, a) -> -> 

^e^a)r{T,a) = ^i?(„)(6)(Q;(e)(T, B)) — — ^e(6)s((/)(T, ,^(t, ct)) , rs(T, ,^(t, ct)) ) , 

'^(a)(T,^) = E / ciVi/C[„),(a,ai,T|a(e),r>,r,] ^ e"^ (t, ^i) 

(^, ^l) + 



Tbs n 



+ j dVi i^(a)(b)(o', ai, r|a(e), r, 0, 7r";,)(r, ai) + 



Since there is no canonical choice of an origin in the 3-diffeomorphism group manifold, the choice 
of a 3-coordinate system is done in two steps: i) by adding the gauge fixings ^''(r, (j) — a'' «i 0, 

declaring that this 3-coordinate system {(X^} on S,- is conventionally chosen as origin of the 3- 
coordinatc systems; ii) by choosing a parametrization of the reduced cotriads, and therefore of the 
3-metric associated with the chosen 3-coordinate system, only in terms of the 3 functions ^(r, a), 
ra{T, a). 
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^5'rs(T, ^) = XI ^e(a)r(T, a) ^6(^)5 (t, ct) = 

a 

ae(r,a) ar(r,a) 



E i\a)u Sa)v] (0(t, e(T, a)), r,(T, ^T, a)) ). (1.4) 



Here are arbitrary rotation matrices and ^e^ay are reduced cotriads depending only on the 
three functions and r^. The kernels K,'[a)ui ^(a){b) ^"^^ ^\a) '^'^^ known explicitly: they 
are the solution of quasi-linear partial differential equations determined by the canonicity of 
the point transformation and by the fact that the rotation constraints do not have vanishing 
Poisson brackets with the super-momentum constraints [3]. Once we have found the general 
solution of the quasi-linear equations for the kernels and we have done a definite choice of 
gauge, the six rotation and super-momentum constraints give further restrictions on the 
kernels, which amount to restrict the general solution to a particular one. 

The family of 3-orthogonal gauges (all having a diagonal 3-metric, ^Qrs = ^e{a)r ^^{a)s = 
for r s), which is parametrized by the last gauge variable 7r<^(r, a), is determined 
by the gauge fixings ^^(r, a) — a'^ ^ and by the following parametrization of the reduced 
cotriads 



'9rs = E 'ka)r 'ka)s = ^rs 6^ (1.5) 
a 

The numerical constants -far, satisfying Es lav = 0, E« laullu = Es lulav = Kv-\, 
define a one-parameter family of quasi-Shanmugadhasan canonical bases (1.3). In the gauge 
ai(a)(T, fa 0, ^'■(r, B)-(t'' ^ 0, we have ^e(„)r(T, a) ^^a)r{T, a). 

The physical deterministically predictable degrees of freedom of the gravitational field 
are the non-local DO's (their expression in terms of the original variables is not known) 
ra(r, a), 7ra(r, a), a = 1,2. In general they are not Bergmann observables, i.e. coordinate 
independent quantities, being non-tensorial and coordinate-dependent. 

Even if we do not know the solution (f) = vr^, rs, tts] of the Lichnerowicz equation, 
the class of Hamiltonian gauges defined by the gauge fixing x(t, a) = 7r0(T, a) ~ has the 
special property that the DO's ra{T, a), 7ra(r, a) remain canonical also at the level of Dirac 
brackets, so that these gauges can be named radiation gauges. When in a radiation gauge 
there is no other residual gauge freedom, we can express every tensor over the space-time 
only in terms of the two pairs of canonically conjugate DO's. 
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This allows for the first time to arrive at a completely fixed Hamiltonian gauge of the 
3-orthogonal type ^, which when restricted to the solutions of Einstein's equations (i.e. 
on-shell) is equivalent to a well defined choice of 4-coordinates for space-time. 

It is evident that the Hamiltonian gauge variables of canonical gravity carry an in- 
formation about observers in space-time, so that they are not inessential variables like in 
electromagnetism and Yang-Mills theory but take into account the fact that in general rela- 
tivity global inertial reference frames do not exist ^. The separation between gauge variables 
and DO's is an extra piece of (non-local) information [13-15], which has to be added to the 
equivalence principle, asserting the local impossibility to distinguish gravitational from iner- 
tial effects, to visualize which of the local forces acting on test matter are generalized inertial 
(or fictitious) forces depending on the Hamiltonian gauge variables and which are genuine 
gravitational forces depending on the DO's, which are absent in Newtonian gravity Both 
types of forces have a different appearance in different 4-coordinate systems. 

In every 4-coordinate system (on-shell completely fixed Hamiltonian gauge) 

i) the genuine tidal gravitational forces appearing in the geodesic deviation equation will 
be well defined gauge-dependent functionals only of the DO's associated to that gauge, so 
that DO's can be considered generalized non-local tidal degrees of freedom ; 

ii) the geodesies will have a different geometrical form which again is functionally depen- 
dent on the DO's in that gauge; 

iii) the description of the relative 3-acceleration of a free particle in free fall given in 
the local rest frame of an observer will generated various terms identifiable with the general 
relativistic extension of the non-relativistic inertial accelerations and again these terms will 
depend on both the DO's and the Hamiltonian gauge variables 

Therefore the Hamiltonian gauge variables, which change value from a gauge to another 
one, describe the change in the appearance of both the physical and apparent gravitational 
forces going (on-shell) from a coordinate system to another one. 



^Namely with ^grs{T-,ff) diagonal and with ^grr{T,a) = fr{ra{T,B)) after the solution of the 
Lichnerowicz equation. The 3-orthogonal class of gauges seems to be the nearest one to the 
physical laboratories on the Earth. Let us remember that the standards of length and time are 
coordinate units and not Bergmann observables [21]. 

^The equivalence principle only allows the existence of local inertial frames along time-like 
geodesies describing the world-line of a scalar test particle in free fall. 

^°When dynamical matter will be introduced, this Hamiltonian procedure will lead to distinguish 
among action-at-a-distance, gravitational and inertial effects. 

^^Note that the coordinate-dependent definition of gravito-magnetism as the effects induced by 
'^gT-r is a pure inertial effect, because it is determined by the shift gauge variables. 
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In this paper we shall start from tetrad gravity in the preferred completely fixed 3- 
orthogonal gauge 7r<^(r, a) ~ and we shall define a background- independent Hamiltonian 
linearization of the theory after this gauge fixation, which, on-shell, corresponds to a definite 
choice of a 4-coordinate system. 

The standard linearization [22-24] of Einstein's equations selects those space-times which 
admit nearly Lorentzian 4-coordinate systems and whose 4-metric is well approximated 
by the splitting ^g^j^^ = "^rj^v + '^h^^, I'^^/ivl << 1 {weak field approximation) with ''r/^j, 
the background Minkowski 4-metric in Cartesian coordinates. Then it is assumed that 
da ^h^y = j^OC^h/j,^), where for the length L one can take the reduced wavelength X/27i of the 
resulting gravitational waves Namely, one assumes the existence of solutions of Einstein's 
equations admitting nearly Cartesian 4-coordinates and split them in a background and a 
perturbation, with a residual gauge freedom in the choice of 4-coordinates. Then, after the 
restriction to the family of harmonic gauges, one replaces the residual gauge freedom on the 
harmonic 4-coordinates with the gauge freedom of a manifestly covariant spin-2 theory over 
Minkowski space-time with Cartesian coordinates. Cravitational waves are usually analyzed 
in the special TT harmonic gauge, a special case of the Lorentz gauges of the spin-2 the- 
ory. Therefore, there is a discontinuity in the conceptual interpretation and one gets that a 
special relativistic theory with its associated absolute space-time chrono-geometric structure 
and theory of measurement replaces Einstein space-time with its problematic concerning the 
identification of point-events (general covariance and the Hole Argument) and with a the- 
ory of measurement in curved space-times still to be developed (we have only an axiomatic 
theory employing test non-dynamical objects) [13-15]. 

Instead, by working in the preferred completely fixed Hamiltonian gauge with well de- 
termined 3-coordinates on the WSW hyper-surfaces and well determined lapse and shift 
functions (so that the reconstruction of Einstein space-time and of its 4-coordinates can be 
done by using the embedding of WSW hyper-surfaces) , we have the possibility of making the 
linearization of the theory on the WSW hyper-surfaces by approximating the Lichnerow- 
icz equation and the Hamiltonian, namely the weak ADM energy, with their linear and 
quadratic parts in the DO's respectively, without never introducing a background. In this 
way we may evaluate the linearized conformal factor of the 3-metric and the linearized lapse 
and shift functions and to obtain a linearized Einstein space-time with 3-orthogonal coordi- 



^^The other gauge fixings are (r, a) 0, Q;(a)(r, a) ^ 0, ^^{t, a) — a'^ and Eqs.(1.5). 

^^With the flat Minkowski background, the background Riemann tensor vanishes and the back- 
ground radius of curvature TZ is infinite {JZ~'^ is of the order of the Riemann tensor). If A is 
the ampUtude of the gravitational wave, the weak field approximation is vahd if .4 << 1. In the 
more general short-wave approximation [22] the background is a vacuum Einstein space-time with 
4-metric '^g\Il^ and typical background radius of curvature TZ. The sphtting '^g^y = ^gj^'^ + ^hfj,,^ 
is done in steady 4-coordinatcs, where, if A is the amplitude of the perturbation, one has: 1) 
'^hjj,i, = 0{A); 2) da'^gjuP = 0{TZ~^); 3) da^^h^y = 0((^)~^). The short-wave approximation is 
valid if ^ << 1 and « Tl. 
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nates on the WSW hyper- surfaces (with a hnearized embedding into the space-time), namely 
a post-Minkowskian Einstein space-time hnearization of a Christodoulou-Klainermann space- 
time. It turns out that this gauge is not a harmonic gauge (and, as a consequence, not a 
TT harmonic gauge), but rather it is a radiation gauge for the gravitational field without 
any residual gauge freedom. The DO's ra{T,a) replace the two polarizations e+ and 
of the TT harmonic gauge. Moreover, besides the absence of background, there is also no 
post-Newtonian expansion: our approximate solution of Einstein equations in the chosen 
4-coordinate system describes a linearized Einstein manifold with Lorentz signature near to 
the Minkowski space-time. It can be called a post-Minkowskian space-time with a linearized 
gravitational field dynamically modifying special relativity. 

In the special 3-orthogonal gauge with 7r0(r, a) = 0, the canonical variables for the 
gravitational field are ra{T,a), TTa{T,a), and (f>{T,a) = e^^"^'^^/^ = 0[ra(T, a), 7ra(T, a)] is the 
conformal factor solution of the Lichnerowicz equation in this gauge. Since the linearization 
consists in taking only the terms quadratic in the DO's inside the weak ADM energy and 
only linear terms in the Lichnerowicz equation, we will do the following assumptions: 

A) We assume |ra(r, a)| << 1 on each WSW hyper-surface and \dura{T,a)\ ~ jO{ra), 
IdudvTaiT, a)\ ~ ■^O(ra), where L is a big enough characteristic length interpretable as the 
reduced wavelength X/2ti of the resulting gravitational waves. Since the conjugate momenta 
7ia{T, f?) have the dimensions of ""^^3°" , i.e. of |; with k = y^;<j) we assume |7ra(r, a)\ ~ jO{ra), 
\duna{T,a)\ ~ j20{ra), \dudyTia{r,a)\ ~ j^O{ra). Therefore, rs(r, a) and Ha{T,a) are slowly 
varying over the length L (for ra,7ra — >• we get the void space-times of Ref. [3]). From 
Eq.(D7) we get that the Riemann tensor of our space-time is of order ^('^s) — J2 ^('^«) ~ 
7?."^, where IZ is the mean radius of curvature. Therefore we get that the requirements of the 
weak field approximation are satisfied: i) A — 0{ra), if A is the amplitude of the gravitational 
wave; ii) ^ = 0(ra), namely « TZ. 

B) We also assume q{T,a) ~ 0{ra), duq{j,(j) ~ j^Oir-g), dudvq{T,a) ~ -^O(ra), so that 
we get 0(r, a) = e'^'^'^''^)/^ 1 + ^q{T, a) + 0(r?) for the conformal factor. The Lichnerowicz 
equation becomes a linear partial differential equation for g(r, a). The linearization is done 
by systematically discarding terms of order 0(r?) in the Lichnerowicz equation and 0(r|) 
in the weak ADM energy. 

With these positions we can, for the first time, solve all the constraints of tetrad gravity, 
Lichnerowicz equation included, and find explicitly the kernels in Eq.(1.4). We can check that 
our gauge is not a harmonic gauge. Notwithstanding this fact, the Hamilton equations imply 
the wave equation for the DO's ra(T, a). We can check explicitly that the linearized Einstein 

equations in this 4-coordinate system are satisfied. Therefore, for the first time we get a 
definition of gravitational wave in a linearized post-Minkoskian Einstein space-time without 
introducing any background and independently from the post-Newtonian approximation. 

After a comparison with other approaches, in which it is emphasized the coordinate- 
dependent nature of effects like gravito-magnetism, we make a study of the time-like 
geodesies and of the geodesic deviation equation. The study of null geodesies allows to 
show the post-Minkowskian modification of the light-cone in each point of space-time. Then 
we analyze the embedding of the WSW hyper-surfaces of our gauge in the post-Minkowskian 
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space-time and the associated congruence of time-like observers. Moreover, anticipating the 
introduction of matter to be done in a future paper, the comparison of our gauge with the 
post-Newtonian approach shows that it is the equation determining the bulk part n(T, a) of 
the lapse function which tends to the Poisson equation for the Newton potential for c — > oo. 

In Section II we solve the linearized Lichnerowicz equation and the rotation and super- 
momentum constraints in our gauge. In Section III we determine the quadratic part of 
the weak ADM energy in terms of the DO's and we evaluate the lapse and shift functions 
of our gauge. Then we find the linearized 4-metric and the tetrads of our gauge and we 
study the Landau- Lifschiz energy-momentum pseudo-tensor in our 4-coordinates. Section 
IV contains a comparison of our results with the Lichnerowicz- York conformal approach, 
with the standard harmonic gauges (and their associated gravito-electric-magnetic analogy) 
and with the post-Newtonian approximation. The Hamilton equations for the DO's and 
their solution, included the verification that Einstein's equations are satisfied, are presented 
in Section V. In Section VI we study the time-like geodesies and the geodesic deviation 
equation of our space-time. Section VII contains the determination of the embedding of 
the WSW hyper- surfaces of our gauge into the space-time and the study of the associated 
congruence of time-like observers and of the null geodesies. Final remarks and the future 
perspectives for the introduction of matter are drawn in the Conclusions. 

In Appendix A we reproduce the results about tetrad gravity [3] needed for this paper. In 
Appendix B we give the determination of some kernels connected with the Shanmugadhasan 
transformation. Appendix C contains the material on the Fourier transform on the WSW 
hyper-surfaces of our gauge needed for the study of the Hamilton equations. Finally in 
Appendix D there is the linearized form of the relevant 3-tensors on the WSW hyper-surfaces 
of our gauge and of the relevant 4-tensors over our space-time. In this Appendix there is 
the explicit verification that Einstein's equations are satisfied. 
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II. LINEARIZATION IN THE 3-ORTHOGONAL GAUGE WITH 7r<^(r,CT) = AND 

SOLUTION OF THE CONSTRAINTS. 



In this Section we study the hnearization of tetrad gravity and the solution of its con- 
straints in the completely fixed 3-orthogonal gauge with 7r^(T, a) — defined in the Intro- 
duction. 



A. The Super-Hamiltonian Constraint. 

Let us first consider the linearization of the Lichnerowicz equation, see Eq.(Al), which is 
the super-hamiltonian constraint of Eq.(A2), interpreted as an equation for the conformal 
factor of the 3-metric. 

Since the super-hamiltonian constraint of Eq.(A2) becomes 7Yfl(T, ct) = Ti^( ~ 

8A0(t, ct) + ■:^Y.aulau9lra(r,a)^ + ■^0(rl) 0, where A — d'^ is the fiat Laplacian 

[A^;^ = —6^{a)], only the term (—A -|- in Eq.(Al) gives a contribution of order 

0{ra). As a consequence the linearized Lichnerowicz equation for (f) — e*/^ = l + ^q + 0(r?) 
becomes 

Ag(r, a) = ^ E7.na^.(r, a) + -lo(r^), (2.1) 
whose solution vanishing at spatial infinity is 

Q{r,a) = --^T.Tau j d'^^ f'f^'] + 0{rl). (2.2) 



B. The Rotation and Super-Momentum Constraints. 

After having solved the super-hamiltonian constraint we have to solve the six rotation 
and super-momentum constraints. As said in the Introduction, this is equivalent [3] to 
find the kernels appearing in the linearization of Eqs.(1.4). After putting equal to zero 
the Abelianized rotation and super-momentum constraints in Eqs.(1.4), the old momenta 

^"^(a) ^) '^—^^(a)i^^ ^)|3-o,7r^=o ^rc givcu by Eqs.(A4), whose linearization is 

W^,)iT,a)^V3Y.Tas [ dVi^[„),(a,ai,T|0,r,]((/)-V7sS5^6.-^7r,)(T,ai) = 

sa 

= ciVi/C[:;;(a,a07r,(r,aO + Ao(r^) = 

= ^3^755 / (iVi[5^5(„)55^(a,ai) + 7^^']^/(a, ai)]7ra(r, ai) + 



^^|3_o means in the family of 3-orthogonal gauges. 
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with 



= - ^ (2.3) 



Therefore, the hnearized kernels /C|°j^'s are determined by the linearized kernels G|°j'^"'s. 

As a consequence, the linearization implies that the partial differential equations (A5) 
for the kernels G[")'s have to be restricted to ra(T, a) — na{T, a) — (zero curvature hmit), 

namely the hnearized kernels G(°)'^"'s are the same as for void space-times [3]. Their general 

solutions G^°J^^'"(a,aiys will determine the /C|°j^'s [Eq.(2.3)] and the F^a)ib)^ [^^^ Appendix 
B]. Moreover, to satisfy simultaneously the rotation and super-momentum constraints , 
the /C(a)^'s must satisfy the linearized version of Eqs.(A6). These linearized equations, 
which restrict the general solution, will be given in Eqs.(2.8) and are equivalent to the 
statement that the old momenta ^T^'^a) satisfy the linearized form of the six rotation and 
super-momentum constraints if we have 

a.%)(T,a) = 0. (2.4) 

Eqs.(2.4) also correspond to the linearization of the three Einstein equations associated with 
the super-momentum constraints of ADM metric gravity. 

The zero curvature hmit of Eqs.(A5) implies that the G|°)^"'s are determined by the 
following linear partial differential equations 

1) s = a homogeneous equations : 

aG|J;^(a,aO^^:r(^,^i)_ 



dal dal 



+ — ^Vr— = 0, a = 1,2, 3; 



daf dal 



^^ri,r2 7^ 6, n 7^ e^a){ri){b) = -^{a)r2^{ri){r2){b)i ^{a){r2){b) = ^{a)ri^{ri){r2){b)^ ^{u){ri){b) 

~^{u)r2^{ri){r2){b)i ^(«)(r-2)(6) = <^(u)n ^(n)(r-2)(6) • 
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dajdaf dafdal da\dal ' 

2)s^as^r,r^a inhomogeneous equations : 

Each set of homogeneous equations, considered as equations for functions of a, is of the 
form d2U{i){a) + diU{2){a) = dzU{2){a) + i92M(3)('?) = f^iM(3)(o') + dzU{i){a) = 0. This is a 
system of three hnear partial differential equations for the three unknown functions it(j) {a) of 
elliptic type, since the determinant of its characteristic matrix [25] is 2^i^2^z 7^ 0. Moreover 
it is integrable, since = /(r)(o"'') with arbitrary f(r) are solutions of the system. We 

do not know whether they exhaust all the possible solutions. Therefore, = 
(7[), with /t(°)"'^ arbitrary functions, are solutions of the homogeneous equations. 

As a consequence, if ai), r 7^ a, is a particular solution of each set of inhomoge- 

neous equations (2.5), then the general solution is ai) = (5|°|'""((7, (Xi) + (7|°j'""((T, ai), 

r ^ a, with the 5'(°)'^"'s arbitrary homogeneous solutions (again with gf|^j'^"(a, cr^), if this is 
the most general solution of the associated homogeneous equations). 

Therefore the general solution of Eqs.(2.5) for the kernels G^^-I^^(a, ai) can be written in 
the following form 

G[:r(a, a,) = (a, a,) + (1 - 5[„))[G|:j™(a, a,) + g^-{a, a,)], (2.6) 

with arbitrary 's and g^a{^^^- Then, Eq(2.3) gives the following expression for the 
kernels /C[^];;'s 

/ oh({ ia.ai) (jG({ (a,ai) og)( ((J,(7i) x 

-fe 'Vy ' +(^-SU)l ' + ' ]). (2.7) 

The solutions of Eqs.(2.5) for the G|°j™'s are restricted by the requirement that the 
^(a)u'^ of Eqs.(2.7) satisfy the zero curvature limit of Eqs. (A6), which in the 3-orthogonal 



17 



gauges become 



^i) 



(o)r 



0, Qy^b, 



= 0. 



(2.8) 



The first set of Eqs.(2.8) becomes the following set of three linear partial differential 
equations to get the Sf^aj^'s of Eq.(2.6) with a 7^ 6 in terms of the G^^j^^^'s 



( (o)au/-> -> \ (o)bu/-> -> \\ 

-( j=m„,K<T). 



(2.9) 



For each pair a ^ b, this is a system of three elliptic linear partial differential equations for 
the /ab's. Each choice of the ^'(^^""'s, a 7^ 6, which gives a solution of this system, implies 

that the associated kernels /C(°)^'s satisfy the rotation constraints. 



Having found a solution for the (?|^j'""'s, a 7^ 6, the second set of Eqs.(2.8) becomes the 



following set of equations for the /i(a)""'s of Eq.(2.6) in terms of the G^^^^'s and 5'(°j™'s 



{o)ruy 



(o)r«5 



= E^r- [(^(a)<^(a)«(^^(o','?l) 

(o)r 
(a) 



(2.10) 



By using the hnearization of the Green function of the covariant divergence given in 
Eqs.(A8) and (A9), we get {f^^^ are solutions of the homogeneous equation) 



= /(T)(^>^l) 



d^<^2 J2 ^isKa)C'{a, ^2) E [Sla)^{a)J%ff2, ^i) 

(a) r 2 



9G\J) ((T2,(7i) 9^;^; ((72, (Jl) 



+ 



(2.11) 



Again this is a system of elliptic linear partial differential equations for the /i(a)"'"s with 
fixed a. 



^^At zeroth order {(f){T,a) = 1) in the 3-orthogonal gauges the spin connection vanishes, 
3w^(„)(r, a)=0 + 0{ra), see Eq.(A7). 
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After having found the solutions of Eqs. (2.5), (2.9) and (2.11) for the kernels G(^)™'s, 



'(a) 

/C|°j^'s and Pj^^jl^^^^s, every remaining arbitrariness will be fixed by the boundary conditions 



at spatial infinity for the momenta ttJ^'^(t, a), given in Eqs. (1.2). The final solutions are 
equivalent not only to the solution of the rotation and diffeomorphisms constraints, but also 
to their Abelianization in the 3-orthogonal gauges with Q;(a)(T, a) = 0. 

C. A Solution of the Rotation and Super-Momentum Constraints. 

We have found the following particular solution G^^^^^ = G^^J^^^ of Eqs. (2.5) with g^^j"'^ — 
^(a)"" ~ vanishing for |a| — > oo (we use S{a^, a{ — > — oo) = for finite cr''), 

G'S^V.^O^O, a = 1,2, 3, 

Gg^^(a, a,) = Gg^^(a, B,) = G^^\B, B,) = G^^\B, a,) = 



1 
2 



G[°f^(a,aO = G^$;\a,a,) = G[f'{a,a,) = G^^;'{a,a,) = 

= I £l dwl S'{B, wlalal) = ^ e{al, a') 5{a\ al) 5{a', af), 

G[°p(a, a,) = Gg"(a, a,) = G^f ^(a, a,) = Gg^^(a, a,) = 

= -2^(^i'^ j^'lf^i,^- ) ^ , 

G^^\a, a,) = G^^;\a, a,) = ^ £l dwl dwf S^a, wla^wl) = 
= "2 ^(^^1' ) — ^1^1' ^ )' 

G'[2)^^(a,ai) = G[3j^^((T,ai) = -\ J_'^dwl J^^dwl6^{a,alwjwf) = 

1 g(5(aSa|) 2 2N 3 3x 
= ~2 ^^'^1'^ J^^(f^i>cr j. 



or 
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2 
d 



0(T\ J —oo J —oo 

[1 - 5(a)] [1- 2(5™ + U] 



d 



(2.12) 



where 9{x) is the step function [-j^9{x) — 5{x)]. This imphes the following expression for 



the kernels /C 



(o)r , 
(a)u 



{a)i 



with 



^i)!:(^'^0=0, a = 1,2, 3, 



^i)l^(o',o^i) 



^(o)2 



(1)3 l^'^^lJ 



= 7^,\a, ^i) = ^ 5^(^, <rlwlaf) = 

1 a^^oj) 2 r/ 3 3x 

= ~2 aa| ^v^^i'^ ^"^'^ ''^1'' ^ 

= (a^ II '^^^^ £ 



1 a 



(2)2 



2 7- 



^fl(^l ^^i^^l^ A(^3 3^ 
"2 <^ ) 5((7 , (7i j 



1 92 



2 (9(7?)2 



-00 J —00 



-00 ^ —00 

32 »„i 



^fl(^l ^l^fl(^2 2^^^5(oVl) 



^(0)1 ^ X 19 



2 5((T ,(Ti)6^((Ti,(T ) = 

o /a n2 / ^^1/ c^w?53(a,WiVi2w?), 
2 (c/fJij ^-00 i-00 
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^(o)i/^ ^ V Id 



l) J —oo J—oo 



^i) = ^ £ £ 
1 a 



2 {dalf 



^3)3 ^1) 



"2 '^('^ ''^i) j = 

1 92 



2 (9(7?)2 



-00 J—oo 

dwld^{a,alwlal) = 



1 d rl 



2 9(7i J -00 



1 ^2 



2 {dalf 



dwl / dwl5''^{a,a\wlw\), 

-00 J—oo 



or 



— /"^^ dw[ r dwtS'{a,w[wta1^^^'^^) = 

X ) J—oo J—oo 



{da 



(2.13) 



The kernels 's are given in Appendix B. 
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D. The Old Cotriad Momenta from the Solution of the Constraints. 



Eqs.(2.3) imply that the cotriad momenta, solution of both the linearized rotation and 
super-momentum constraints, have the following expression in terms of the DO momenta 

^7r(i)('r, a) = ^7ai7ra(r, a). 



^'^) = [7al J ciVi 



^(ai,cr )c)(cr ,crij7ra(r, cri) + 



2 7ra(T, (7ij 



v5 
2 



r 1"°° 



+ 



+ 7a2 / da^ — h 



- 7a3 / da^ / da^ 



1 r ^■2d''7r-aiT,alala^) 



2 



-la2 j d ai9{ai,a )d{a ,cti) — na{T,ai) 



+ 



+ 



roo roo 



+ 



+ 7a3 / dai 



'7f(2)(r,a) = 7r(i)(r,a), 

^7i'(2)('r> ^) = \/3 E 7a27ra(T, ct). 



3-3 



V3. 



j3^ <92_5(aVi)^^^2 ^2^^^^3 ^3^ 

{daiy 



7r(2)(r, ^) = E [7.1 / d'a.'^^^^-^eiat, a^ial, <7^)7r,(r, ^0 + 



.1 .l^^^(^^^ ^1)^(^3 ,3> 

+ 7a3 



+ 7.2 / d'a,S{a\ al)^^^^^e{al a')7r,{T, B,) + 



d5{(7\ erf 



T^a{r, ai) 



V3 



T.[Ta. I da! I da! 



.d^7 ra{T,a^a!a! 
{da^Y 



+ 
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iiO ^ 








POO 










+ 7a3 / 


1 


3-1 






3-3 . 




^(3) 


[r, 


^) 


= 




qx2 






3-3 / 




^(3) 




^) 






3-3 
^(3) 


[r, 


^) 




7a37r 



or 



a au 



Clearly Eqs.(2.4) are satisfied. 

The solution (2.12) for the G(°j™'s is such that the momenta ^7r[„)(T, a) of Eq.(2.14) tend 
to zero for \B\ — > oo, as required by Eqs.(1.2), if the momenta TTair, a) satisfy the restrictions 



+00 

ci(7''7ra(T,a) = 0. (2.15) 

00 



For instance these restrictions are satisfied if TTa{T,a) = q^iq^q^z with ^g^f^^ ^(t*->oo 
[i-ij-ik cyclic] in a direction-independent way. 

We have not succeeded in finding a solution without these restrictions. Eqs.(2.15) can 
be thought of as 6 additional constraints defined on 2-dimensional surfaces. As we shall 
see, the consistency of these restrictions with the final Hamilton equations will impose the 
following restrictions on the DO's ra(T, B) 



+00 

da'r-a{T,a)^0. (2.16) 

—00 



Therefore Eqs.(2.15) and (2.16) are 6 pairs of second class constraints and we could 
think to find a Shanmugadhasan canonical transformation from the DO's ra(r, cr), 7ra(r, c?) 
to a new basis in which 6 pairs of conjugate variables vanish due to Eqs.(2.15),(2.16) and 
the physics is concentrated in the remaining pairs. However, we shall not look for such a 
transformation, because it is highly non-trivial due to the fact that these constraints are 
defined only on 2-dimensional surfaces. We shall go on to work with the DO's ra{T,a), 
TTair, a) even if this will imply formal complications. 

In conclusion, we have been able to solve all the Hamiltonian constraints of tetrad gravity 
on the linearized WSW hyper-surfaces of our gauge. 
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III. THE WEAK ADM ENERGY, THE LAPSE AND SHIFT FUNCTIONS AND 

THE 4-METRIC. 



After the solution of all the constraints, in this Section we determine the weak ADM 
energy, the lapse and shift functions and the 4-metric of our linearized theory. Since we 
need the weak ADM energy, namely the Hamiltonian in the rest-frame instant form of 
tetrad gravity, see Eq.(All), we need the following results of Appendix A: Eqs.(AlO), where 
there is the expression of the weak and strong ADM Poincare' generators, and Eqs.(A12), 
where there is the expression of the weak ADM Poincare' generators in our gauge in the 
canonical basis (1.3). 

In Eqs.(All) we put: i) (j) — 1 + + 0(r?) with g(r, a) given by the solution (2.2) of 
the linearized Lichnerowicz equation; ii) the expression (2.14) for the cotriad momenta. In 
this way we get the form EADM[ra, '^a] of the weak ADM energy only in terms of the DO's 
of our completely fixed gauge. 



A. The ADM Energy and the Lapse and Shift Funcions. 



The Hamiltonian linearization of tetrad gravity in our completely fixed gauge is de- 
fined by approximating the weak ADM energy EADMl^a, 7i"a] with the quadratic functional of 
ra(r, a) and of na{T, a) contained in it 

Eadm = -^Padm,r = —3— / d^'^ Z]^s('r > ^) 
2AnG 



+ 



+ 



EE<^ra)7»"7to / d^ai£a2r^^f{ai,a2)7ra{T,ai)7r-h{T,a2) + 

ah 



+ ^EE7..75. / d'aSa,d'a2[Y.Tl:^:{a,a,)Tl^:{a,a2) 

ah '^s « 



+ 
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E ^^^-^ E7 
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d 



ab 



4:n\a — ai\ 



+ 0(4). 



(3.1) 



See Eq.(C13) of Appendix C for the Fourier transform of the ADM energy. 



The determination of the bulk lapse function n{T, a) is done by solving the integral 
equation (A13) written in our completely fixed gauge. Therefore we must evaluate ^"^^1^^'^^^ 

and From the linearized version of the super-hamiltonian constraint given before 

Eq.(2.1) we get 



SHr{t, ai) 



— e- 



5(l){T,a) "27:0 
while from Eqs.(Al2) and (2.1) we get 



(3.2) 



e 



6(j){T, a) 



dVi [25ir-?(T, ai) 7=Y^ardirra(r, ai)\dir5^(a, <Ti) + 0(r?) = 

V3 a 



— e 



2TrG 
+ O(rI). 



Ag(r, a) - ^j:Tardya{r, a)] + 0{rl) 



(3.3) 



Then the integral equation (A13) for the lapse function becomes the following partial dif- 
ferential equation 

3 

j ciVin(r,ai)Ai5=^(a,ai) = -e^An(T, a) = + O(r^), 



27rG 



n(T,a) = + O(r^). 



(3.4) 



The determination of the hulk shift functions uses this result and the linearized form of 
Eq.(Al4) (see Eq.(C14) of Appendix C for the Fourier transform), 



Y: I dV2[C<^(„)„5^(^l,^2)+7^1)l"'(ai,a2) 



l]7a^;7ra(T, a2)G'j^j"''(ai, u) + 0(r?) = 

y ;>(t" ra°' poo roD 

\ Y i^-^i^uv + Sav)] dal dal dad da", 
aV,(r,a,V,Vr"'") 
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/(T^ ("cr" /-oo /-oo 

u^r '°° °° '^i '^1 



In conclusion in our completely fixed gauge the lapse and shift functions are 

N{t, a) = -e + n(r, a) = + 0{rl), 
Nr{r, a) = Urir, a) = -e '^Qrvir, a). 



(3.5) 



(3.6) 



Since the shift functions are of order 0(ra), we get the following results : 

i) the 4-coordinates associated to our Hamiltonian gauge are not synchronous, so that we 
are using non-time-orthogonal reference frames and wc cannot use Einstein's convention for 
the synchronization of clocks [26]; a non-standard definition of simultaneity of distant time- 
like observers is needed, consistent with the Hamiltonian description based on the Cauchy 
simultaneity WSW space-like hyper-surfaces E,-; 

ii) there may be coordinate-dependent gravito-magnetic effects; 

iii) the velocity of light becomes direction- dependent (sec Rcf. [26]) : if -u* is a unit 3- vector 
with respect to the 3-metric ^jrs — — e(^gfrs+ ^Ig" ), i-e. ^7rsM'^M* — 1, [see Appendix A of Ref. 

[2], after Eq.(A5)], the light velocity in direction is w{u^) — y„^^'" ^ — l—u^nr + 0(r?) 




B. The Linearized 4-Metric. 

After the solution of all the constraints and the determination of the lapse and shift 
functions, the /^-metnc of our linearized space-time with S^'^'^'^^ -adapted coordinates in the 
3-orthogonal gauge with 7r0(r, a) = becomes (we write it in the form of a perturbation of 
the Minkowski metric in Cartesian coordinates only to visualize the deviations from special 
relativity of this background-independent post-Minkowskian space-time in the 4-coordinates 
associated with our preferred 3-orthogonal gauge) 

^9ab{t-, a) = '^riAB + *^/iAs(T, a), 



4 



hrAT,a) = + O{rl), 



av 

av 



hrr\T, Cr) = -enr{T, a) ^ - e — AGEM,Z-0,rV^^ ^) = 3-^ 2^7, 

[ [l-2(5„„ + 5a.)] 

ua,u^a 

/ct" /-(t" i-oo rOD 

dal / dal / dal^ / da.. 



f)2 ( u a kj^u,a\ 
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2^[l-2(5„, + 5,,)] 

/c'^ ra^ POD roo 

da{ j da'^ / / da. 

-oo J —oo Ja^ J 



) + 



+ 0{ri), 
-2e 



(li^^^) + larTair, B) 6rs + 0{rl) 



.J^i^^) = -6^7(r, a) + O(r^) = -6 06(r, a) + O(r^) = -e [1 + 3 g(r, a)] + O(r^) 



— e 



1 - 



d^ai 



47r|a — ai| 



Y^(r,a) = e + 0(r^), 

Y^(T,a) = -e5-n,(T,a) + 0(r^), 



^r^(r, a) = -e ( 



— e 



47r|CT — CTil 



(3.7) 



Therefore, the 3- metric on the hnearized WSW hyper- surf aces of our gauge is only con- 
formal to the Euclidean 3- metric, namely the linearized S^^'^^^ of our gauge are conformally 
flat. 

The linearized cotriads, triads and adapted cotetrads and tetrads become [see Eqs.(A7) 
and (A16)] 



2 j 



,3^ dlur-a{T,Bi] 



47r|(j — (?i 

^e(„)(T, a) = [l - ^ I] (7arra(T, a) + 



+ 



e(r, a) = J7(^, f?) = 1 + 3g(T, a) + O(r^) 
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47r|a — (7i| 



+ 0(r^). 
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grs{r,a)^-e grs(r,a)^ e(^a)r(r,a) e^a)s(r,a) ^ 



^l^lau J a ar 



47r a — (Ji 



(3.8) 



4 IT' ^ :t\ lA 



UE^o)(r, a) = l^(r, a) = -6(1; -5^^n,(r, a)) + O(r^), 



(E)%)(r,a) = (0;=^e[,)(r,a)) 



= (0;<^[a)[l-^E(7..r.(r,a) + 



x(o) 



(T,a) = /^(r,a) = -e(l;0) + O(r^), 



x(a) 



47r|(T — (Til 



]) + 0(rD. 



(3.9) 



Since in our gauge we have 99(a) (r, a) = 0, the tetrads and cotetrads "^E^^ and "^E^^ of 
Eqs.(l.l) coincide with those of Eqs.(3.9). 



Sec Appendix D for the associated Christoffel symbols, spin connection, field strength 



and Riemann tensor of and for the 4-tensors of M^. 



C. The Energy-Momentum Landau-Lifschitz Pseudo- Tensor. 

Usually the energy of a gravitational wave on the Minkowski background is evaluated 
[22,24] as a mean value over various wave-lengths of the coordinate-dependent Landau- 
Lifschitz energy obtained from the Landau-Lifschitz symmetric pseudo-tensor [27]. 

The Landau-Lifschitz pseudo-tensor (L)t'^'^ = {L)t^^ ■, which contains no second derivatives 
of the metric and gives meaningful results only in an asymptotically fiat Cartesian coordinate 
system, was found trying to rewrite the consequence ^ ^T'^^ = of the Bianchi identities 
applied to Einstein's equations in the form + = 0. Starting from Einstein's 

equations '^T'^" ^ 'Sg^^^^^ ~ l^fl'''^^-^)) o^i® rewrites them as 

(L)^'''' + \l^g ^T^" = dph^^^p^ , with 
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Then one gets 

iovrG z 

+ ^(2 V - V^)(2 ^5 V/j - V^7/3)9p V'^S. Y"]. (3.11) 

However, as noted in Ref. [28] and emphasized in Ref. [29] , if the energy flux carried away 
to infinity by gravitational radiation in a given asymptotically fiat space-time is deduced 
from the Landau- Lifschitz pseudo-tensor, then the result is reliable only on an appropriately 
chosen flat background metric and then a connection between Bondi energy at null infin- 
ity and ADM energy at spatial infinity can be established [29]. Since in our approach to 
Hamiltonian linearization of Einstein's equation we have no background metric, the Landau- 
Lifschitz pseudo-tensor is not a useful quantity. 

If we write the Landau-Lifschitz in S^'^'^'^^-adapted coordinates, (i)t^^(r, a) and then 
we choose our special 3-orthogonal coordinates, then the Landau-Lifschitz 4-momentum 
Pll = / d^'^ {L)t^'^{T, a) has to be contrasted with the weak ADM energy Eadm, namely the 
Hamiltonian in the rest-frame instant form of gravity, and with the vanishing weak ADM 
momentum Padm ~ evaluated in those coordinates. 

/^From an explicit calculation done by using Eqs.(3.7) it turns out that does not 
agree with Eadm and Padm [for its expression see Eq.(5.11)] 



5 5 
- E 9jAdjfi - E drfidiUi - o E QiAdifi + 

i,j i i 

+ ^ '^{dsTlrf + ^ E dsUrdrUs - E dsUsdrTlr + 
r,s r,s r,s 



+ 0{rl) ^ Eadm, 
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5 5 1 

1 7 

~ 2 ^ + drfrdrfr " 3 «9iAajnr + 

4 

+ o X] drAdiUi + ^ diAdrUi - ^ difrdiUr - 

i i i 

+ o{tI) + p;^^, 



A(r,a) = x/3^7.n/ciV^ 
o.^ J 47r cr 



2 

= ^Il7aira(r,a). (3.12) 

Therefore, when in a future paper we will add matter to tetrad gravity, we will have to 
devise a method independent from the Landau-Lifschitz pseudo-tensor to identify the energy 
of the matter and its variation due to the emission of gravitational waves. 
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IV. COMPARISON WITH THE LICHNEROWICZ-YORK CONFORMAL 
APPROACH, WITH THE STANDARD LINEARIZED THEORY IN HARMONIC 
GAUGES AND WITH THE POST-NEWTONIAN APPROXIMATION. 



Since in the literature there are many coordinate- dependent definitions of gravito- 
magnetism, which are a source of ambiguities, in this Section wc shall review same of them 
and we will rephrase them in the language of our linearized post-Minkowskian space-time. 
Moreover we will show that our completely fixed Hamiltonian gauge does not belong to 
the family of the harmonic gauges used in the standard background-dependent lineariza- 
tion of Einstein's equations. Finally we will make some comments on the post-Newtonian 
approximation of our space-time. 



A. Comparison with the Lichnerowicz-York Conformal Approach. 



To establish the connection with the Lichnerowicz-York conformal approach [18,19,16] 
we need the extrinsic curvature of the WSW hyper-surfaces of our gauge. From Eq.(A15), we 
get that the extrinsic curvature and the ADM momentum of the WSW hyper-surfaces of our 
gauge in the hnearized theory are (see Eq.(C15) of Appendix C for the Fourier transform) 

^ u 

= eVS^^ [2dr slasT^ air, B) + 
^ a 

+ [1 - Srs]Ylaw[i - "^iSrw + Sgyj)] 
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^ rs ^ rs 



(4.1) 



since the trace, proportional to York's time T = — e^|^ '^K, turns out to be 

'Kir, a) = 5- 'Krsir, a) + 0{rl) = + O(r^). (4.2) 

This means that the linearized WSW hyper- surfaces are constant mean extrinsic cur- 
vature (CMC) surfaces and satisfy the maximal slicing condition. Moreover, the ADM 
super-momentum constraints are satisfied. 

Therefore, York's traceless distorsion tensor ^A^^ = ^K^g ~ \ 'drs "^K = 'Krs + O('^i)) see 
Ref. [16], Chapter 4.10 and Appendix C of Ref. [30], coincides with the extrinsic curvature 
in the hnearized theory. After a conformal rescahng 'grs — 0^ ^Qrs we get 'A^^ — ^A^^ + 
'A^£ = 0^° 'A''" ~ 'K^^ + 0{rl), where 'A^j. and ^A^* are the transverse traceless (TT) and 
longitudinal (L) components respectively. The longitudinal component satisfies ^AY\s = 

^A'-'\/=^ {AlWY = AW"^ + l{W'\sy^ + ^K'sW' , where W''{T,a) is York's gravito -magnetic 
vector potential. In the hnearized theory it satisfies A{61 + = dg^K^^ + 0(r?) = 

+ 0(r?), whose solution is 

W^{T,a)^0 + Oirl), 



3 Ars (^LW^^ ^ 

+ O(r^). (4.3) 



3 3 



Therefore, York's transverse traceless physical degrees of freedom of momentum type, 
namely 'Ai^j, = 3^'i _3^rs _ ^]^rsj^Q^^2^^ linear functionals of the vTa's of our linearized 
theory. The coordinate-type physical degrees of freedom in 'g^s depend upon the r^'s 



B. Comparison with the Standard Linearized Theory in Harmonic Gauges and its 
Associated Gravito-Electro-Magnetic Analogy. 

Let us remark that our special 3-orthogonal gauge is not a member of the family of har- 
monic gauges, because, in coordinates adapted to our gauge, the condition Ba = 
becomes 



^^All the quantities evaluated with the rescaled metric having det \'g\ = 1 will be denoted with a 
^^In the 3-orthogonal gauge we have 'grs = ^gfT^ , see Eq.(Al7). 
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dA {^Y^)] {r, a) = [d^ (v^ - l) - E ris] (r, a) + 0{rl) = + O(r^), 

s 

dA {\f^^f^)] {r, a) = [dr rir - dr (-^ Y^barra){r, a) + 
+ ^ Et../ ^) + 0{rl) ^ + Oirl). (4.4) 

Therefore, on the solutions of Einstein's equations [15] our radiation gauge identifies 
a 4-coordinate system different from those compatible with the standard linearized theory 
coupled to matter. 

This theory, in the post-Newtonian approximation, allows to make the gravito- electro- 
magnetic (GEM) re-formulation of Einstein's theory (see Refs. [16,31]), emphasizing the spin 
1 aspects of this spin 2 theory over the Minkowski background in Cartesian coordinates. 

By putting '^g^^ = '^t]^^ + h^^ with /i^^ = h^^ + Yv^^h {h = V V, V^V*" = 
'^T]^^ — h'^", h^^ = '^rj^°' '^Tj^l^ hoiis) Einstein's equations in presence of matter become □ h^i, = 
^Tj"^ dadph^y= — Y[.T^jj when the Lorentz gauge condition d^h^^ = 9^ (^^/^^^?^^) = 0; 
viz. the conditions for harmonic coordinates (which include the post-Newtonian ones), are 
imposed. Even if in this approximation there is no back-reaction of the gravitational force 
on the sources, it allows to make an analogy with electromagnetism [23]. h^i, replaces the 
4-potential with the coordinate transformation (changing the Christoffel symbols but 
not the curvature tensor) h^^, = h^y — d^h^ — dy -\- '^rj^^ da corresponding to the gauge 
transformation = -|- 9^ A. Maxwell equations □ A^^ — dfj da A" = — ^ correspond to 
the linearized Einstein equations □ h^i,-\- '^r]^y da dp If^ — d^ da h^^ — du da h"^ = —^T^i,. The 
Lorentz gauge dnA^ = (with residual gauge freedom given by functions A's satisfying □ A = 
0) with equations OA^^ corresponds to the harmonic gauge (harmonic coordinates) 

dt, h^'^ — with □ h^i, — —\ (with residual freedom in the choice of coordinates having 
the 6^'s satisfying nb'^ — 0; in Cartesian coordinates the equations of motion are decoupled 
one from the other). 

The retarded solution h^u{x°, x) = ^ J d^x' '^^"^^ ^'^^J^ ^'^ ^ allows to define: 

i) the gravito- electric or Newton potential hoo = — e 4 ^gem determined by p = T°° /(? 
as the effective gravitational charge density; 

ii) the gravito-magnetic vector potential (the shift functions) of Ref. [16,31] hoi — 
e ^ Aqemi determined by = T^'/c as the effective gravitational current density; 

iii) to disregard at the lowest order the 3- metric hij = 0(c~^). 

As shown in Ref. [31], in the harmonic gauge the gravito- electro-magnetic 4-potential 
i^GEM, Agem) satisfies the Lorentz gauge condition 

do^GEM + d- {^Agem) = 0. (4.5) 

The factor ^ derives from the fact that the effective gravito-magnetic charge {Qb = source) 
is twice the gravito-electric charge {Qe — ^smirce)- The gravito- electro-magnetic fields are 
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Egem = -d^GEM - doi^AoEM), Bqem = 9 X Aqem = CQ.GEM- in the harmonic gauge 
(harmonic coordinates) they satisfy Maxwell equations 

^ ^ - 1 - 

d ■ Egem = 4:TtGp, d ■ {-Bgem) = 0, 

d X Egem = -^^(^-Bgem), d x {^Bgem) = doEGEM + -^J' (4-6) 

and one can consider as electro-magnetic-like gauge transformations the residual coordinate 
freedom. 

In this approximation the GEM fields are determined only by the matter and vanish in 
absence of matter, viz. they are independent from the physical degrees of freedom (like Va, 
TTa) of the gravitational field. 

If we consider solutions of the homogeneous equations Ohfj^i, = 0, d^h'^'' — 0, they are 
used to describe gravitational waves h^^ = a^,^e*'^"^" (a^^jy = const., kfj_k^ = 0, a^^.k'^ = 0) as 
the only independent degrees of freedom of the gravitational field. The residual freedom of 
coordinate transformations in the harmonic gauge may be used to get a^^ = a^jy — c^kj^ — 
Cv kfi + '^rj^i, Ca with a^^, = a'^j^^ = and a'^^ = a^* = through the choice = —i c'^ e*'^"^" 
(solution of □ ft'* = 0). In this completely fixed harmonic gauge (with non-diagonal 3-metric), 
named TT gauge, only two independent degrees of freedom (polarization states) survive: 
they are the counterpart of r^, tts of our completely fixed 3-orthogonal radiation gauge with 
7r0(T, a) — 0. If k'^ — (0; 0, ^, ^), the only non-zero a^^ of the transverse gravitational wave 

are a^^, d^y = — a^^, d^y = a'y^, and we get ds'^ = e {dx°Y — (1 + fxx) dx^ — 2 f^y dxdy — (1 — 

fxx) dy"^ - dz"^]^ with e+ = f^^ = a'^^ cos {z - x°) + v?), Cx = fxy = cos (^ {z - x°) + ^). 
In this gauge the lapse and shift functions n, Ui are all zero and there is no GEM potential 
{^gem, Agem) coming from the transverse gravitational wave. 

Let us remark that in the electro-magnetic case on Minkowski space-time with Cartesian 

— * — * — * — * 

coordinates the radiation gauge Ag — d-A — 0, with the transverse fields A±, E± as the only 
physical degrees of freedom (DO's), is a particular case (name it gauge T) of the Lorentz 

gauge d'^ Afj_ = 0, obtainable with a A = —ice^^'^ with k'^ = from □ A = 0. Instead, the 4- 
coordinate system on the post-Minkowskian space-times connected with our radiation gauge 
for the gravitational waves is different from the 4-coordinate system which is reinterpreted 
as the TT harmonic gauge of the spin-2 theory on the background Minkowski space-time 
with Cartesian coordinates. 

Instead in our non-harmonic radiation 3-orthogonal gauge with 7r^(r, a) = 0, even in 
absence of matter the shift functions do not vanish and we can define 

i) a vanishing gravito- electric potential f?) = + 0(r?); 

ii) a gravito-magnetic vector potential Aq'^Ii{t, a) = {— e ^^^(t, a)} determined by the 
shift functions; 

iii) a gravito-electric field Eq'^Ij{t, a) — —dr [\ A%'^1^[t, a)] — —e^dr nr{r, a) also de- 
termined by the shift functions. 

iv) a gravito-magnetic field BQ'^l[{T,a) = d x AQ'^lj{T,a) = cQ,Q'^l[{T,a), with 
^^^'^(t, a) — {f^GBM.r- — <^~^ ^ruv |^ } the gravito-maguctic precession angular velocity 
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(connected with the prcccssional effects of Lense-Thirring or dragging of the inertial sys- 
tems) again determined by the shift functions. 

Since we are not in a harmonic gauge the Lorentz condition is not satisfied d- (| ^^^m) 7^ 
0. The analogy with Maxwell equations is partially lost in our completely fixed (radiation) 
gauge: while the equations deriving from the existence of the potential d ■ (l^^^'^f) = 
and d x Eq^^^ = —dr{\ B^gem) hold, the other two are not satisfied: d ■ e'^q'em 0, 9 x 
{\Bq'^Ij) 7^ QtEq'^Ij. This shows how the gravito-electric- magnetic analogy is coordinate 
dependent. 

Finally, the coordinate transformation (a passive 4-diffeomorphism) from the E^.^'^^^- 
adapted 4-coordinates r, a of our completely fixed gauge to the 4-coordinates = a;^(r, a) 
of the standard harmonic TT gauge and the relation between the DO's ra(r, a) and the TT 
polarizations e^, can be obtained as solution of the following system of partial differential 
equations: 



4 / \ dXy 4 , 4 

QixAx) = - — gAB[r, a), g^^ = 



feO \ 

— ee+ — ecx 

— ecx +ee+ 
V -e / 



(4.7) 



with '^gAB given by Eqs.(3.7). 

If we write x^{T,a) = fij.{r,a) + gf,{T,a) with = 0(1) and = 0{ra) and if we 
require the asymptotic behaviour /^(r, a) -^\ff\-,oo X{oo)^t{0) +Uf,{PADM) t + er ^{u{Padm)) o'' 
with w'^(p) = p^/^/ep^, then we get that 

i) the four equations with "^goo = e, '^goi = can be solved to give ^ in terms of 

ii) the system of partial differential equations^g^ia = ^^^23 = 0, ^gzz — — e and ^gu — —^g22 
have to be used for the determination of the ^-dependence of the four functions ^^^^(r, a); 

iii) the use of the solution of point ii) inside point i) allows to find ^ and then, by 
integration, the functions a;^(T, a); 

iv) finally the equations e+ = —e'^gn and = —^*gi2 give the two polarizations e+, Cx 
in terms of the Dirac observables ra(r, a). 



C. Connection with the Post-Newtonian Approximation. 

The standard post-Newtonian approximation [32] is applied to the gravitational field 
created by an isolated compact object, like the Earth, described by an energy-momentum 
tensor T^''' [T°° = 0{c^), T°' = 0{c), T'^ = 0(c°)]. It is a weak field {^g = ^ri + % near 
zone approximation giving corrections to Newton gravity. 

Following Ref. [33], if = is the Newton potential of the compact object (of mass M 
and radial dimensions Rg < R), the post- Newtonian (PN) expansion of the 4- metric ^gf^u is a 
series in the dimensionless parameter C ^ ps ^ (T is a characteristic time of variation of 
the source). The IPN approximation keeps the following terms in "^gnu- ^goo = e + ^ + 

9oi — 9ij ~ ~l~ "3^"- 
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In Ref. [32] the PN approximation was given in a PN 4- coordinate system (adiabatic 
approximation of a Robert son- Walker line element appropriate to a homogeneous isotropic 
cosmological model containing the compact object), defined as a local, quasi- Cartesian 4- 
coordinate system at rest respect to the universe rest frame and in a standard PN gauge 
where the 4-metric has the form (we use the notation of Yi^J^^'^) -adapted coordinates) 

V. = 6e-^^ = 6[i-^y + o(c-^)], 

V. = -^K + 0(c-^), 

V. = -e^Qrs = -eSrs (1 + + 0(c-^), (4.8) 

with U = V m. the standard PN gauge. Here V is the Newton potential generated by the 
compact object. 

Now our 4-coordinate system, in our completely fixed Hamiltonian gauge, is by definition 
in the rest frame of the universe (rest-frame instant form). In presence of matter, the 
linearization condition (as we will show in a future paper on tetrad gravity plus a perfect 
fluid, where it will be shown that both the lapse and shift functions depend on matter so 
that n 7^ 0)) requires T^^ = 0{ra) and we have 

^g^^ = e [(-e + nf - ^Qrsrfn'] = e(l - 2en[matter\) + 0(r?), 

^Qrr = -^^grsrf = -enr[r a, matter] 0(r?), 

Vs = -^^Qrs = -e(^rs (1 + KVa, matter]) 0(r?), 

0=[det^^]Vi2 = (l + ^) + O(r^). (4.9) 

Therefore, in our gauge the Newton potential is the lapse function^ V = n[matter] (and 
the equation determining the lapse tends to the Poisson equation), the gravito -magnetic 
potentials are the shift functions, Vj. — Ac'^ nr[r a, matter], and the Lichnerowicz equation 
for the conformal factor 4> amounts to a determination of U, which may not coincide with 
the Newton potential like it happens in the standard PN gauge. By solving the linearized 
equations for the DO's (see next Section), Dra = Za[matter], i.e. by considering their 



^^In Refs. [34,33] there are various definitions of 4-coordinate systems. The standard PN gauge is 
said to correspond to the conditions dg ^Qts ~\dr ^Qss = 0{c~^), dg ^Qrs — \ {^gss — '^9tt) = 0(c~^). 
Instead the algebraic spatial isotropy condition of Ref. [33], —^gff^grs = ^rs + 0{c~^), contains 
both the harmonic and the standard post-Newtonian gauges. Finally a ADM Hamiltonian gauge, 
needed to include 2.5PN gravitational radiation reaction, is = 0, dg ^grs ~ \ (^r ^9ss = 0. 
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Lienard-Wiechert retarded solution without incoming free radiation, we can recover a form 
like the one of Eqs.(4.8), which does not depend explicitly on the DO's 

Let us now consider the relation of our linearized radiation gauge with Ref. [35], where 
the Galileo generally covariant formulation of Newtonian gravity, depending on 27 fields, 
was obtained as a limit c — oo on the ADM action of metric gravity. 

The final action of Ref. [35] with general Galileo covariance depends on the 26 fields Ao, 
ao, Ar, ar, ^Qrs, Jrs, Prs (after having put the field = 1 by rescaling the absolute time) . 
There are 18 first class constraints and 8 pairs of second class ones. It turns out that cko, Ar, 
ar, three components of ^Qrs, one component of the momentum conjugate to ^Qrs, the trace 

and the longitudinal $^ parts of $rs in its TT decomposition, and the longitudinal 7^ part 
of are Hamiltonian gauge variables, while A^, (the Newton potential) and the remaining 
components of ^Qrs, 7rs, Prs are determined, together with their conjugate momenta, by 
the second class constraints. There are no propagating dynamical degrees of freedom. The 
gauge variables describe the inertial forces in arbitrary accelerated non-Galilean reference 
frames. 

To make a comparison with the results of Sections 6 and 7 of Ref. [35] the starting point 
is the following parametrization of the 4-metric of Eq.(3.7) (we show only the 26 terms which 
appear in the Newtonian action) 

9rr-e[i-^ + ^ + 0{c )] = 

^e{l-2en) + 0{rl), 
v. = e [A + ^ + 0{c-')] ^-enr + O(r^), 

grs = -e grs = -e^^, + — + — + 0{c ) = 
= -eSrs{l + kr) + 0{rl), 



n[matter] ^ + + 0{c ^,r?)]. 



c2 c4 
1 



nJra, matter] ~ —Aj. -ar + 0(c , r=). 



Let us also remark that strictly speaking to reach Newton gravity one usually performs a double 
limit: i) the zero-curvature limit ra,7ra 0, eliminating the genuine degrees of freedom of the 

gravitational field, which do not exist in Newtonian gravity; ii) the 1/c expansion. However the 
zero curvature limit is not consistent with Einstein's equations with matter and must be replaced 
(in a linearized theory) with the restriction of the DO's r^, tTq to the Lienard-Wiechert solution. 
Then the 1/c expansion will kill the curvature for c — ^ 00. 
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(1 + kr[r-a, matter]) = + + ^ + 0{c-',rl). (4.10) 
Instead the PN approximation (4.8) implies 

Ao = eV, Ar = 0, ar = -4K-, ^Qrs = Srs, Irs = U 6rs, (4-11) 

which are consistent with the gauge freedom of Ref. [35] (it is a possible gauge of the general 
Galileo covariant description of Newtonian gravity) . 
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V. THE HAMILTON EQUATIONS AND THEIR SOLUTION. 



In this Section we shall study the Hamilton equations in the preferred 3-orthogonal 
gauge and we will verify that their solution produces a hnearized solution of Einstein's 
equations. Since, as said in the Introduction, the admissible WSW hyper-surfaces must 
admit a generalized Fourier transform, wc shall use it on our linearized WSW hyper- surfaces 
(where it coincides with the ordinary Fourier transform on at the lowest order) to study 
the Hamilton equations. 

As said in the Introduction, in the rest frame instant form of dynamics on WSW hyper- 
surfaces the Dirac Hamiltonian is H]j = Eadm + Mt) • Padm- However we can add the gauge 
fixings JYdm ~ [i-^- vanishing of the weak ADM boost generators of Eq.(Al2)] to 
the three first class constraints Padm ~ with the consequence A(r) = 0. This means that 
we have to eliminate the internal 3-center of mass of the universe (it is put in the centroid 
■^{oo)i'^) ^^^^ spatial infinity used as origin of the 3-coordinates on WSW hyper-surfaces), 
i.e. 3 pairs of global degrees of freedom among the DO's ra(T, a), 7ra(r, a), reducing them to 
canonical variables relative to the internal 3-center of mass. 

Like in the case of Klein-Gordon and electro-magnetic fields and like for every isolated 
system treated in the rest frame instant form [5], we should find the canonical transformation 
from the canonical variables ra(r, a), 7ra(r, a) to a canonical basis containing the canonical 
internal 3-center of mass Qadm: weakly equal to the M0ller internal 3-center of energy 
^ADM — ~j'ADM/-^ADM) Padm{^ O) as the conjugate momentum, and the internal relative 
variables R^, Ha- Usually [36,37], like in the Klein-Gordon case, we should start with a 
naive 3-center of mass X^dm^ conjugate to Padmi we should find the canonical variables 
i?^, relative to it and then we should use the Gartenhaus- Schwartz transformation to find 
[36] the canonical variables with respect to Q^adm- Then we should add the gauge fixings 
Q\dm ~ P-ADM — ~JYdm/^adm ~ 0, which put the internal 3-centers in the centroid x^^^ 
origin of the 3-coordinates on the WSW hyper- surfaces and implies Hu = Eadm, to the 
rest-frame condition Padm ~ we should go to Dirac brackets with the result Ra = R^, 
Ha = H^. Presumably also the final relative variables have to satisfy Eqs.(C2). 

However, since also the determination of the naive 3-center of mass Xadm of gravity is 
not trivial, we will fix this final gauge freedom simply by putting equal to zero the Dirac 
multipliers, A(t) = as a pre-gauge fixing condition compatible with Jadm ~ 0- 

We can now study the Hamilton equations for the independent canonical degrees of 
freedom ra{T,a), 7ia{T,a) of the gravitational field, only restricted by Eqs.(2.15), in our 
completely fixed 3-orthogonal gauge generated by the Hamiltonian Eadm of Eq.(3.1). These 
Hamilton equations replace the equations n/i^i/ = 0, d^.h^'^ = of the standard linearized 
theory in the harmonic gauge coordinates. 

A. The First Half of Hamilton Equations. 

The first half of the Hamilton equations associated with the Hamiltonian (3.1) yield the 
following expression for the velocities drVair, a) 
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drVair, a) = {ra{r, a), Eadm} = ^^^7ra(r, a) + 



c3 



u,r,b,(a) 



l2nG 

+ 



c3 

u,r,s,b (a) 

+ E ('^(a)^'^(6)« - '^(a)«'5(fe)t,) / d?(Tirl-°f{a,ai) j (iV27^^''(<Ti,a2)]7r6(r,a2) + 

^(a)(6) 

+ \0{rl). (5.1) 

By using the expression of T^°y given in the last two hnes of Eqs.(2.13) it can be 

checked that the two terms hnear in the T's vanish due to a factor 5(a)u(l — 5(a)«). After 
some calculations it turns out that the two terms bilinear in the T's give the same result, 
so that the final expression for the velocities is 

9r'^a(T,^) = ^^^7rs(r,CT) + ^^E^arTfes E - '^{Kr + 5^rW- - '^{Ks + ^vs)] 

brs uv,uj^v 

/■0-" r<^Y /-oo /-oo 3*717 (t n'^n'"n^^'^'^'\ 1 

L /: L k ^ i'"^-'- 

In obtaining this result we have made integrations by parts justified by the asymptotic 
vanishing of 7ra(r, c?) and we have used dx6{x,y) = —dy6{x,y), dxS{x,y) — —dy5{x,y) as- 
sumed valid on linearized conformally flat WSW hyper-surfaces. 

To invert these equations, to get the momenta in terms of the velocities, we shall assume 
the validity of the Fourier transform on the linearized WSW CMC-hyper-surfaces. They 
are conformal to R^, i.e. they are Euclidean plus corrections of order 0{ra), irrelevant 
when acting on the functions ra(r, a), 7ra(r, cr); for the scalar product we have ^cjrsk^h^ = 
k ■ h + 0{ra)- As a consequence functions of order 0{ra) are considered as functions over 
at the lowest order of approximation. 

By using the Fourier transform defined in Appendix C, Eqs.(5.2) become Eqs.(C4) 

drf-a{r, ^) = E A,dk)h{r, k), (5.3) 

6 

with the matrix Ag^j^{k) and its inverse given in Eqs.(C5) and (C6), respectively. Let us 

— # 

remark that, notwithstanding ^ab(/c) diverges for k'^ — > 0, Eqs.(5.3) are well defined if 
7fa(T, ^) — (A;^/c^/c^)^+^ /a(^), 6 > 0. This condition is stronger of the requirement (C2), 
i.e. of the Fourier transform of Eqs.(2.15). 

As a consequence we get [e^^ = — e^^, ejg = 1] 



40 



7ra(T, k) = ^ A_5^(^)9^fa(r, k), 



7ra(T, ^)^Ylj d^CTlGabi^ - ai)drn{T, CTi), 



^3 



A-^{k)e'''<^-^'^ = 



0^5/1 TT: 1 N n / ^ ON n / ^ — 5T7r + 



. ^ \- ^OA IVOA W(w3-^l) 

2 ..Ye/''^'''''"^' " " hdw^,ndw^,ndwir 



The final result of the inversion is 



8nG J (47r)3|(T — ai \ \ai — |ct2 — CT3I 

r d^drrajr, as) 

ka<73i)2(a<7|)2(a<7i)2 + 

2 ,,ftc-/''^'''''"^' " " ^(a^smiHa^ 



(5.4) 



(5.5) 



Eqs.(5.5) satisfies the condition (2.15) automatically. 



B. The Second Half of Hamilton Equations. 

Let us now study the second half of Hamilton equations associated with the Hamiltonian 
(3.1). After some calculations we get 

drTTair, a) = {7rs(r, a), Eadm} = 

^ d-i\, {T. a) (9^5 (r, a) 



The Fourier transform of Eqs.(5.6) is given in Eqs. (C7) and (C8), with the remarkable 
result d-r7rair,k)= - A'^iik)fair,k). 

This implies da^ Vair, a) = 0, i.e. Eq.(2.16), as the simplest way to get consistency 
between Eqs. (5.6) and (2.15). 

Let us remark that, as shown in Eqs. (D9) of Appendix D, our Hamilton equations imply 
the satisfaction of the remaining Einstein equations ^Rrs{T, a) = 0. 

Eqs. (5.4) and (5.6) imply 
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dlvair, a) = '^^^^drTiair, a) + 

+ ^E7ar7-6s E [l-'2{5ur + Svr)]['i--'2{Sus + Svs)] 

brs uv,u^v 

/cr" /"Cr" POO POO 

da\ dal da^J da',{d',d',drn){T,alala'^^'^) + 
-oo J—oo J a'^ J a\ 
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Ur-a{T,a) = [dl-Y.dl]r-a{T,a)^ 

r 

^iv (A - \ [ ^^^^ d^n{r,ai) 
^'^ ' 47r|^ - ^i| {dal^idaff 

+ J EE E(-2(^6c + 476t7a - 76t7c«;)7ar-76.(l - Suv)[i - 2{5ur + Syr)] 

/tr" pa^ POO POO 

da^ / dal / da^ / ^^^2 
-oo J—oo Jct" JfTj* 

/ " ^3 / , , + 7i<-'('^aJ- (5.7) 



As shown in Appendix C, Eq.(ClO), all the terms in the second member of Eq.(5.7) 
cancel, so that we find that the DO's ra(T, a) satisfy the wave equation. Actually, as shown 
in Eq.(ClO), the Fourier transform of the equation of motion is 



□ r^(r,a) = 0. (5.8) 

Therefore, in our radiation gauge we get the wave equation □ra(T, ct)=0 for the DO's. 
Notwithstanding the presence of gravito-magnetism, which, as said, should imply the 
anisotropy of light propagation, we get an isotropic propagation of gravitational waves in 
the radiation gauge. 

The complicated form (3.1), (C13), of the weak ADM energy shows that in our gauge we 
do not have the conventional description of a massless spin two particle over Minkowski 
space-time like in the TT harmonic gauge. This fact is connected with the lack of a 
coordinate-independent notion of gravitational energy density, which has no counterpart 
in the theory of massless spin two particles in Minkowski space-time. Instead we will see in 
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Eqs.(5.11) and (5.12) that the weak ADM 3-momentum and angular momentum do admit 
a standard particle interpretation. 

The solutions of the Hamilton equations for the DO's ra(r, a) and 7ra(T, a) are given in 
Eq.(Cll), (C12) 



(5.9) 



with the functions Ca{k) vanishing as {k k k^Y for k'^ — > 0. 

C. Special Solutions for the Background-Independent Gravitational Waves. 

— * 

Special solutions are the following plane waves whose 3-momentum h cannot lie in any 
coordinate plane [h ^ 0; (1, 0, 0); (0, 1, 0); (0, 0, 1); (1, 1, 0); (1, 0, 1); (0, 1, 1)] 

rair,a) = e^^'^'^-l'^l"), C-a = const., 

f-a{r,k)^{2nfC,e-'\^^\^S\k-h), 
n.ir, a) = -i \h\ ^ Azifi) C,e^^^-'-\^'\^\ 

b 

n-a{T,k)^-i{2nf\h\Y,A-{{h)C-,e-'\'^\^S'{k-h), 7fs(T, 0) = 0, 

b 

(5.10) 

which verify Eqs. (2.15) and (2.16). Here tt^ has been evaluated with Eq.(5.5) and Eqs.(C4) 
and (C7) are valid also at k — 0. 

Hoever these solutions do not satisfy the rest frame conditions P^dm ~ 0) which restrict 
the solutions to globally outgoing or ingoing wave packets. 

Before looking for these solutions, we shall give the form of the weak ADM charges 
(A12) in the radiation gauge with our solution of the constraints (see Eqs. (3.1) and (C13) 
for Eadm)- The weak ADM 3-momentum and spin of Eqs.(A12) assume the same simple 
form in terms of Fourier transformed quantities as for free masslcss fields due to the fact 
that we are in an instant form of dynamics. This is not true for the ADM boosts, which 
have a complicated form like the ADM energy. 
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For the 3-momentum, due to exact cancellations, we get 
Padm = - d^<^Yl l^r rc(r, a) - 

•' c 

-1-cri-T.Tau I C^Vi ^y.^^'t] -2j2rar Or r,{T, a))] 7T,{t, a) + 
^ au J 47r|(7-(7i| s 

cs 



+ E 4^ E Tam J d'a, + E 9. r,(r, a)) + 



/7T V o /am ^t) I ^ . i ^ ^ 



Ml) V <-< ~ am 



V3 am 47r |(7 - (7i| ^ ^ ^ 

^ - d^a J2 '^c{r, o") <9^ rg(r, a) = 

•' c 



I 



' Y^n,iT,k)k^f,{T,-k), (5.11) 



(27r) 



where we used '7^^°^" = 0, see Eq.(2.13). 

For the angular momentum, due to similar cancellations, we get 
Jadm = d^'^Yl ^cir, o) {a^ ds - a' dr) n{T, a) = 

^ 5 

= - / E Ur, k) ik' ±, - k- A) ,,^r, -£). (5.12) 

Eqs. (5.11) and (5.12) reflect the fact that we are in an instant form of the dynamics. 
Instead for the boosts we get a complicated expression hke for the energy (see Eqs. (2. 13) 
for the kernels) 

a ^ ^ am ■> 47r|(7-(7i| ^ 

+ 2(X: r^drn{r^a)f)] - ^ [2 E + 

I ^ a 

+ 4 E lira n{r, B) ( d'a, E ^0 E 7.m 7r,(T, a,) + 

bu ^ a 
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E ^2) E ^5 ^2) + 

till ma 

E ^6?(^'^2) E rbsn{r,^2))]) - 

SttG J 3 ^ 4 y 47r|cr-cri| 

+ E (7a^; - 7a«) ?^a(T, a)) 
a 



^ an J 47r\a-ai\ ^ ' 



+ Ajd'a,Jd'a, E ^it(^,^i)EM(^'^i)^i?(^''?2)E7c-.7r^^^^^ 

u,m,s,(a) b ^ 

f . 2 11 

- y E('^r,; - 1)(2? + E7arra)(2 ^r-? + J^iTbv " 76r)5rr6)(r, a) + 

+ 0{rl). (5.13) 
The weak ADM 3-momentum of the solutions (5.9), (Cll), has the form 

Padm = 2 E / 7^ A-cs'ik) \k\ K Ciik) cm « 0, (5.14) 

be ^"^^l 

so that its vanishing is a condition on the Fourier coefficients Ca{k), which, for instance, 
cannot be satisfied by the plane waves (5.10). 

A class of solutions of the Hamilton equations with the correct vanishing behaviour at 
spatial infinity, which, besides Eqs.(C2), also satisfies the rest frame condition Padm = is 

obtained by taking coefficients C-aik) = Ca ^^^tIt^ e"^' (C^ = const.) in Eq.(5.9) 

|fc| 

?f\-C- [ {kik2hf p r i(\k\T-k-a) 



4Ca 


(1 + 


a 


2_^2) 


(27r)2 9Vi 9V2 [1 + (1^ 


+ rr] [l + {\a 


-r)2] 



d^k {kik2k^f 
(2^ 



■Ka{T, a) = -I / ——- ;^ 6 
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\k\i:^-jik)c. 



g-i(|fe|r-fe-CT) _ gi(|fe|r-fc-CT) 



:47rG 



nr{T, a) = i a/3— ^ 7s^(l - 5„c)[l - 2((5„„ + - 2(5„r + 5cr)\ 



c 

(2^ 



I kj.k^ 
T2T2 



YA;i{k)C-,sin{\k\r)e''-', 



Eadm 


— Yj ^ab CoCfc, 

db 


pr 


= 0, 






fr 

''ADM 


"2 "^AIJM — 



(5.15) 



ab 



The two constants Cj, C2 have to be expressed in terms of the two boundary constants 

M = Eadm, S = \Jadm\ defining the mass and spin of the post-Minkowskian Einstein 



space-time. The solution (5.15) goes hke \a\ 



00, i.e. much faster than the 



^ for |a| 

behaviour (1.2). It should be possible to find solutions saturating Eqs.(1.2), namely such 
that by using Eqs.(3.7) we should get ^Qrs -^r^co (1 + ^1°^ )^rs, = —e^Qrr -^r^oo 0{r~'^). 



By comparison, in presence of compact matter the expected solution should have a post- 



1 - ^, 'g^ 



(1 + ^)S,. 



Newtonian behaviour ^Qrr -^r- 
(see Eq.(19.5) of Ref. [22]). 

Therefore, in absence of matter, the rest-frame condition destroys the transversality prop- 
erty of the TT harmonic gauge plane waves. 
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VI. THE TIME-LIKE GEODESICS, THE GEODESIC DEVIATION EQUATION 
AND THE POLARIZATION OF GRAVITATIONAL WAVES IN THE 

RADIATION GAUGE. 



A. The Geodesic Equation. 

Let us consider a time-like geodesic with affine parameter s, ct^(s) = (^r(s); (t"(s)^ . 
After the hnearization in the radiation gauge the geodesic equation becomes (see Eqs.(C16) 
of Appendix C for Fourier transforms) 

= -a.n.K.))(^) -2 ^rMs))^^- r.,K.))^- — ^ + 

+ 0{rll (6.1) 

where Eq.(D5) has to be used for the Christoffel symbols. 

If we parametrize the geodesic as cr^(s) = + h^s + /^(s|ra, tTq], where a"^ + &^ s is the 
flat geodesic, we get at the lowest order ^ -b'h' ''f;^((T(s)) + 0(r?) -b'h' ''f ;^(a + 

6s) + 0(r?). Therefore the solution to the first of Eqs.(6.1) is 

r(s) = a" + ¥s - /' dsi H ds2^t:Ja + 6s2)&"6" 

^ ^ r ds2^f ;„(a + 6s2)fo"6" with ^f;„(a + &S2) ~ 0(rs). (6.2) 

as JO 

Since the tangent p"^(s) = ^^^^j^ to the geodesic must satisfy cither '^g'^^ PaPb = e (time- 
like geodesic with s as proper time) or '^g'^^ PaPb = (null geodesic) at the initial time, the 
constant b'^ is determined in terms of the 6'"'s and of the initial data. Therefore we will have 
— bj^g-^ + bj^^^ with bj^^^ = 0{ra)- Therefore Eq.(6.2) may be rewritten in the form 

r(s) = + + &[^] s - dsi J^' ds2^tUa + 6(o)S2)6''6" 
P^is) = ^ = &(o) + ^(1) - /J ds^'rUa + b^o)S2) b^b\ (6.3) 



The spatial part of Eqs.(6.1) becomes 

-dru^ia + 6(o)s)(6[„))' - 24f;^,(a + 6(o)s)6[„)6'- - 3f;f,(a + 6(„)s)6'-6* + 0{rl) 



ds^ 



a"(s) = a^ + b^-J^ ds, £ ds2 [(6[„))2 - 2 ;!,(a + 6(o)S2)6[o) 6'' - 'C,(a + 6(<,)S2)6''6' 
P"(^) = = ^" - io - 2'r;^,(a + 6(o)S2)6[„) 6'' - (a + 6(„)S2)6^-6^ 



(6.4) 
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In the time-like case at s = 0, where t(0) — a^, we have the condition 



\s=0 



1 = (6-)2 _ 26^6^71^ (a^,rs (a")) - [S^ - e^hrt{a\ra{a^))b'b% ats^O, (6.5) 
whose solution is 

= + = 6^n,(a, r-a{a)) ± sJ{¥nM, r-a{a-))Y + 1 + [5,* - e4/i,,(a-, r-a{a-))] ¥b\ 



2 v^TO^ 



= r^.(a^ r,(a-)) 6^ ^ J '^"^^ , + ^(r?)- (6-6) 



If we choose b^ — 0, we can take the solution /"^(s) = so that we get r{s) — a'^ + 
>6T=i s + a^. 



With the choice b'^ — 1, — 0, — 0, we get r(s) = s and ^ ^Ja^^^ = — (?t^«(t(s), = 
a + /(s)) = at the lowest order, with the solution /"(s) = — J"^ ds nu{s\ a^) + 0(r?). 
Therefore, with 6"^ = 1, 6^ = 0, a"^ = we get the following expression for the geodesies 

a^{s) = (r(s) = s; a"(s) = - [ dw n„(w, a'')) + 0(r2), 

^ = (1; -M.,a^))+0(r^). (6.7) 
This is the trajectory of a massive test particle in our coordinates. 



B. The Geodesic Deviation Equation. 

If we consider two nearby freely falling particles following two nearby geodesies (6.7) 
such that crj^(O) = (0; 0) and o"^(0) = (0; a") with an infinitesimal a", we get for Ax^(s) = 
a^{s) - a^{s): Ax"(s) = 0, /\x'\s) ^ [5"" - /o^ciw %y V "- If «" = '^"''a is their 



coordinate distance, then their proper distance is A/ = y — e ^S'ab Aa;"^ Ax^ = Ax*" y/^Qrr = 

(1 + iA;,) Ax" + 0{rl) ^ [1 + \kr - Jo dw ^"''^^;""^ ] a + 0(r?) by using Eqs.(3.7) [in a TT 
harmonic gauge one gets A/ ^ [1 + hJ^F] a if r is one of the two polarization directions]. 

As a consequence the connecting vector Ax^(s) satisfies 

This is the geodesic deviation equation along the geodesic o^l{s) of the form (6.7): it 
shows explicitly the action of the tidal forces in our 4-coordinate system. More in general, 
by using the equation 
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ds ds ds 

+ 'r^^{s,a-)Ax'^is) + Oirl), (6.9) 



ds 

the geodesic deviation equation [22,24] takes the form 



ds^ 



+ 2^r^^(.,a")^-^ + 



d^Ax^{s) 
di^ 



= + O(r?), ^ Aa;^(s) = cs + 



+ 2 



(is ds 



1 4x^7" / 

+ Z Ax^(s) = -^^V,,(s, a") Ax^(s) + 0(r?). (6.10) 

ds 

If we take Ax'^(s) = as it happens for the geodesic (6.7), i.e. c — d — 0, the geodesic 
deviation equation becomes 

ds^ ' ds ds ^ ' 

= -^R\sr{s, a^)Ax'{s) + 0{rl). (6.11) 

If we write Aa;'"(s) = A(o)a;''(s) + A(i)x''(s), with A(o)a;''(s) = (o; a") the flat deviation 

and A(i),T'(s) = {q-- !^dw[n^{w,a-) -n^iwM) ~ (O; - /o' c^w ^^^^^U.=o a'') = 0(rs), 
Eq.(6.11) becomes 

/a ^r.x , A ^r. x\ _ A(i)X'-(s) _ 



{A^^x^{s) + A^)x\s)) 



ds^ ^ ' ' ' ' ds^ 

--2 r,,(s,a) A(,)x (s)- 

-^^^,.(s,a")A(„)X^(s) + 0(ri). (6.12) 

The use of Eqs.(D5) and (D7) allows to check that Eq.(6.12) coincides with Eq.(6.8). 
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C. The Tidal Forces Generated by the DO's on Test Particles. 



Eq.(6.8) can be used to see the effect of the tidal forces generated by the two DO's 
ra(r, a) of our radiation gauge, replacing the two polarizations of the TT harmonic gauge, 
on a sphere of test particles surrounding a test particle sitting in the origin (r = 0; = 0) 

at s = and described by the geodesic cr^[s) = (s; (t"(s) = — Jq dwn^{w, 0) ). If we define 
= ar]^ with ^QuvV^V^ = l + 0(ra) and a << 1 (the radius of the sphere), the test particles 
on the sphere at s = will follow the geodesies cr^(s) = (s; cr"(s) = ar]^ — Jq dwn^{s, at]'") ). 
The connecting vectors Ax^(s) = A(o)X^(s) + A(i)a;^(s) = a^{s) — cr^(s) satisfy Eq.(6.8), 
namely 

We have to solve these equations with the shift functions (3.5) evaluated on the solutions 
(5.15), which are not transverse plane waves due to the rest frame conditions (this forces us 
to consider a sphere of test particles) . 

The quasi-Shanmugadhasan canonical basis (1.2) adapted to our gauge is not unique, 
because there is still the freedom associated with the numerical constants jar (their variation 
amounts to a redefinition of the DO's). As a consequence of the conditions stated after 
Eqs.(1.5), the 7ar's niay be written in term of an angle ip with < < 27r 



7lr 
l2r 



'Cosip simp cosip sinip 2 



(- 



cosip sinip cosip sinip 



cos 



(6.14) 



By defining the DO's. i.e. our polarizations, with the convention = 0, we get 71,. = 
^ (1, 1, -2), 72, = ^ (-1, 1, 0). Then Eqs.(3.7) imply 



-e'gSir, B) = 5rs + 'hfXr, Bp 5^^ [1 + A;?)(r, a)] + O(r^), 



k^{T, a) = k^{T, a) = -e— n r, a - - / rfVi — '■ — — — — 

6 2 J 47r I (7 — (Ji I 

Ur, = X [-2ri(r, a) - - j d a, ^—-^ ], 



-e'g?J{T, a) = 5,, + ^/ig)(r, apSrs [1 + k^ir, a)] + 0{rl), 
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kl{T, a) 



-e ^ 3 h ir,a) + -Jd a, ] 



47r|(T — (Til 



(6.15) 



If wc consider our (non transverse) solution given by the wave packet (5.15) and we solve 
Eqs.(6.13) numerically for a sphere of particle at rest around the origin of the 3-coordinates 
on a WSW hyper-surface, we obtain the two 3-dimensional deformation patterns replacing 
the usual 2-dimensional ones for the polarization in the TT harmonic gauge: 

i) in figure 1 there is the deformation pattern for the case Ci 7^ 0, = 0, namely for 
ri(r,a) ^ 0, r2(r,a) = 0; 

ii) in figure 2 there is the deformation pattern for the case Ci — 0, C2 ^ 0, namely for 
ri{T,a) = 0, r2(r,a) ^ 0. 

In the two figures are reported the snapshots at three different times {t = —1, —0.5, 0) of the 
sphere of particles originally at rest (bottom) and the time evolution (from t — —3 to t — 3) 
of the six particles at the intersection of the three axes and the sphere of particle (top), 
whose initial 3-coordinates are (1,0,0) and (—1,0,0) on the x-axis, (0,1,0) and (0,-1,0) 
on the y-axis, (0, 0, 1) and (0, 0, —1) on the z-axis, respectively. Since a particle on a x axis 
will remain on the same axes during the evolution, only the x coordinates are reported for 
the two particles on x axes. The same representation rule has been applied for the particles 
lying on the other axes. 
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FIGURES 
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Particles at t = -1 



Particles at t = -0.5 



Particles at t = . 




D. A Relativistic Harmonic Oscillator as a Resonant Detector. 

Let us remark that in a TT harmonic gauge the basic ideahzation of a test resonant 
detector on the Earth for gravitational waves is a non-relativistic damped harmonic oscillator 
[22,24]. Given two equal masses m located along the r axis at the positions ri(t) and r2{t), 
connected by a massless spring with spring constant k, damping constant v and unstretched 
length lo, one considers the proper extension l{t) — J^'^l^l^ dt y^l + h^^{t) [1 + | h^^{t)\ lo 
in the metric of the gravitational wave and defines ^ — I — lo ^ T2 — ri —lo-\-\ hj!^{r2 — ri) . 
This leads to modify the equation of motion ^ + 2^^ +0;^^ — (u^ — 2k/m, ^ — v/m), 
valid in absence of gravitational radiation, to 

This same equation can be obtained as a consequence of the tidal force at the lowest order 
from the geodesic deviation equation [24] , by considering the center of mass of the detector 
as moving along a geodesic and by identifying ^ with the spatial part of the connecting 
vector Aa;'^(s) with /\x'^{s) = is twice the connecting vector from the center of mass 
to one of the masses) . Note that due to the additional force terms the nearby world- line is 
no more a geodesic. 

Let us try to generalize this resonant detector to special relativity and to sec the effect on 
it of our post-Minkowskian gravitational wave. As shown in Ref. [36], there is a description 
of the relativistic 2-body problem on arbitrary space-like hyper-surfaces E^, leaves of a 
3-1-1 sphtting of Minkowski space-time. If z^{T,a) is the embedding of T,r into Minkowski 
space-time, the two positive energy scalar particles are described by 3-coordinates 77[(t), 
i = 1,2, on E^, such that a;f(r) = z^{T,f]i{T)) and by the conjugate momenta ^^^(t). 
When the hyper-surfaces are the Wigner hyper-planes, orthogonal to the total 4-momentum 
of the 2-body system, we know how to replace the canonical coordinates fji, with: i) 
the canonical internal relativistic center of mass and its conjugate momentum /?_|_ ^ 
(it is vanishing since we are in the rest-frame instant form of dynamics); ii) relative 
coordinates p and tt. Even if we do not know explicitly how to express the old basis r/j. 
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Ki in terms of the new basis due to the non-hnearity of the canonical transformation, we 
are able to evaluate the Hamiltonian for the relative motion in the rest frame. In the free 



case it is = \J ml + 2t{'^ + ^Jm^ + ^ii'^. We can add action-at-a-distance interactions 
either inside the square roots (Droz Vincent - Komar - Todorov (DVKT) models [38]), H — 



nil + ^(/^) + 27?^ + J 1712 + ^(/^) + 27?^, or outside them (hke the Coulomb potential for 



charged particles [39]), H = \Jm\ + 2??^ + ym^ + 2??^ + U{p^). Due to the square roots 
we cannot find explicitly the Lagrangian and the Euler- Lagrange equations for the relative 
motion in the rest frame, except in the special equal mass case. For mi = m2 = m we get 

[36] L — —m \J A: — p' and an interaction of the DVKT type is introduced by the replacement 
M — ^Jm? -\- V{p^). The Euler-Lagrange equations are 



dr 



or 



For V — —^kf? (relativistic harmonic oscillator), the non-relativistic hmit of these equa- 
tions becomes p + o;^ p = 0. A damping term —u p may be added by hand. 

Like in the non-relativistic case, we can adapt this relativistic model, in the case of a 
test 2-body system, to curved space-time: i) the relativistic center of mass moves along 
a geodesic F, with 4-coordinates o"^(s); ii) the relative variable p is the connecting vec- 
tor in St- to a nearby (non-geodesic) world-line Fi. However, since is a Rieman- 
nian 3-manifold, the fiat vector p has to be reinterpreted as the field of tangent vec- 
tors to the 3-geodesic joining F to Fi on each S,- [40]: if 6*^(5, C) is a 3-geodesic with 
affine parameter C, such that F is ^''(s,0) and Fi is ^^(s, 1), then we replace V{p^{t)) 
with y(7^2(s)), where n^{s) = dC^^^^grs{T,e{s,C))^^^ (it is independent from 
C because the tangent vector — is parallel transported along the 3-geodesic). Analo- 
gously we make the replacements p • p 1— > TZi{s) — Iq d( ^^-^^^grs{T,9{sX)) ^ qIqI^^ 

p^ - 7^.(.) = /o^ rfC ^ V.(r, ^(., 0) 

By identifying p''(r) 1— > Aa'^{s), ff{T) ^ ^^^7^, P ^ ^^'^^^■> by adding the force 
terms (damping included) of Eqs.(6.17) in the geodesic deviation equation (6.9) we get 



d-'Ax^s) , r , 4-7^2(5) , 1^ dAx\s) 



+ 



r;.(^) + , ^^ V {n\s)) (^ + - UM) 



+ 



ds 



+ 



d'Ksis) +4^._(^) - ^^^^#TT^'(7e^(5))l Ax^(s) = 0. (6.18) 



ds "--V"/ m2 + T/(7e2(s)) 



21 + 3ri_^{e{s, 0) = with the 3-ChristoffeI symbol of Eqs.(Dl). 
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This is the post-Minkowskian counterpart in our coordinates of the non-relativistic equation 
(6.16) of the TT harmonic gauge. 

For very small |p| we can make the approx- 

imations /c ~ — — |^=o g-rsK'^ [■s)) — — k=o, K-i ~ — — K=o grs[<^ [s)) q^q^ |^=o, 

7^2 ~ ^ agac^'' lc=o ^9rs{c'^{s)) ^ QgQ^'^^ |c=o- this is a pole-dipole approximation of the test 
2-body problem. 
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VII. THE EMBEDDING OF THE WSW HYPER-SURFACES, THE ASSOCIATED 
CONGRUENCE OF TIME-LIKE OBSERVERS AND THE NULL GEODESICS. 



A. The Post-Minkowskian WSW Hyper-surfaces. 

The WSW triads ^e|^'^'^'*'^(T, ct) — >|o!|^oo <^(o)' P^-rallel transported from spatial infinity 
on the hnearized WSW CMC hyper-surface, and the associated cotriads ^e^^^^^^ (r, a) 
are the solution of Eqs.(Al8), whose linearized form is 

L + + ^^f^J (r, a) + 0{rl) = + O(r^), (7.1) 

where Eqs.(Dl) have been used. 
These equations imply 



^\a) ^) = ha)Ar. ^) + ^ ^a) Tan j d a, , 

di fl2)± + ^3 + d2 /^3)± = + 0{rl), 

f(a)±=SU + {S'' + ^)9ia)s: (7.2) 

with gr(a)s arbitrary. The simplest solution is to take g(a)s = 0, so that the WSW triads are 



3 {WSW)r _ rr 



\/3 ^ / 3^ a|^^(r^n 
+ ^ Z^lau a ai — — — — . (7.3) 

2 an J 47rC7-(JiN 



Eqs.(Al9), allow to find the associated WSW adapted tetrads (preferred ADM Eulerian 
observers or asymptotic fixed stars giving a local compass of inertia to be compared with 
local (Fermi- Walker transported or not) gyroscopes) 



4 U'^SW)A 



Once the WSW triads are known, Eqs.(A20) and (A21) give the embedding of the lin- 
earized WSW CMC hyper-surfaces into the linearized space-time 



56 



F'ir, = 'efS'^^^^lT, B)5^a)r(y' + eTn,(T, a), (7.5) 



with the transition coefficients — — g^x^ solution of the hnear partial differential equa- 
tions 

In adapted coordinates Eqs.(7.5) of the embedding of linearized WSW CMC hyper- 
surfaces become 

4wsw){r, ^) = <)(0) + F\t, a), (7.7) 
with (t, c?) = x\oo){^) — ^T- The time-like evolution vector is ^(V5vk)t('^' ^) ~ ( 

{W. 
(a) 



B. Congruences of Timelike Observers. 

In Eqs.(3.9) we have the expression of the E^^'^^^-adapted tetrads and cotetrads. In par- 
ticular there is the expression for the contra- and co- variant normals, 1^{t, a) and Za(t, a), to 
'£(}vsw)^ The associated 4- velocity field defines a (non-rotating, surface forming) congruence 
of time-like observers orthogonal to using the parameter r, labeling the leaves of 

the foliation, as evolution parameter Let us remark that r is not in general the proper 
time of any observer of the congruence. 

As for any congruence, we have the decomposition {Pab — ^Qab — IaIb) 

^Va Ib — Iacib + Pab + ctab + ^ab^ 

— B , 4 — acceleration, 



^^This is the hyper-surface point of view according to Ref. [41]. Instead, the threading point of 
view is a description involving only a rotating congruence of observers (like the one which can be 
built with the 4- velocity field associated to ^(Vsw)r(^' '^)) • since this congruence, being rotating, 
is not surface-forming (non-zero vorticity), in each point we can only divide the tangent space in 
the direction parallel to the 4- velocity and the orthogonal complement (the local rest frame). On 
the other hand, the slicing point of view, originally adopted in ADM canonical gravity, uses two 
congruences: the non-rotating one with the normals to S,- as 4- velocity fields and a second (rotating, 
non-surface-forming) congruence of observers, whose 4-velocity field is the field of time-like unit 
vectors determined by the r derivative of the embeddings identifying the leaves St- (the so-called 
evolution vector field) . Furthermore, it uses the afHne parameter describing the world- lines of this 
second family of observers as Hamiltonian evolution parameter. 
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© = ^Va l^, scalar (volume) rate of expansion, 

oab = ^ {flA Ib + dB I a) + ^ (^Va Ib + ^Vs Z^) - ^ © -Pas, 

rate of shear tensor [with magnitude — -(Tabcf^^), 

^AB — —^BA — ^ABCD ~ 2 ("^^ ~ 2 ^^^"^ ~ ^ ^^-^ ^ ^' 

twistor vorticity tensor, cu^ — e^^'-^^ oubc — vorticity vector. (7.8) 

© is the expansion (it measures the average expansion of the infinitesimally nearby world- 
hnes surrounding a given world-hne in the congruence), uab the shear (it measures how an 
initial sphere in the tangent space to the given world-hne, which is Lie transported along 
l^, will distort toward an ellipsoid with principal axes given by the eigenvectors of g^b 
with rate given by the eigenvalues of g^b) and ujab the twist or vorticity (it measures the 
rotation of the infinitesimally nearby world-lines surrounding the given one); gab and cu^s 
are purely spatial [gabI^ = i^abI^ = 0). Due to the Frobenius theorem, the congruence 
is (locally) hyper-surface orthogonal if and only if cuab — 0. The equation i 9^ Z = | © 
defines a representative length / along the world- line of l^, describing the volume expansion 
(contraction) behaviour of the congruence completely. 

The linearized acceleration of the observers vanishes (i.e. at the lowest level we get 
inertial observers) so that there is no gravito-electric force on test particles in the sense of 
Ref. [41] (geodesic in local rest frame) 

By using Eqs.(D5) we get "^WaIb = -^^^abIc = e'^^AB + 0(r?). The expansion of the 
congruence is 

e = ^Va = dAl^ + Mb = ed-n-ed-n + 0(r?) = + 0(r?). (7.10) 

The shear is not zero so that there is gravito-magnetism on test particles in the sense of 
Ref. [41] (geodesic in local rest frame) 

CTAB = I {aAh + asU) + I CVaIb + 'VbU) = e^T^^B + O(r^), (7.11) 
while the vorticity vanishes (the congruence is surface forming) 

COAB = I {aA Ib - aB I a) + \ {^V a Ib - 'Vs /a) = + 0{rl) . (7.12) 

If Xg^ (r) is the time-hke world-hne, in adapted coordinates, of the observer crossing the 
leave E^^ at we have 
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<(r) = (r;p^„(r)) = 



Xffoi^o) = {ro]^o), A?o(^o) = ^o- (7.13) 

The effective trajectory psoij) is determined by solving the equations p^X^) = <?o + 
Z{wsw){t, PSo{t) ) - Z{wsw){'ro, ^o) with Z(wsw) given by Eq.(7.7). The 4- velocity xj^r) = 

dx^ (r) . ^ 

""^ satisfies 

liir) = i^(r,p,„(r)) = ^t^s ^ , -c < ^^''^^ 

with the observer acceleration given by a4^(r) = ^%;r^ = + 0(r|), a4^(r) /^^^(t) = 0. 

Since Xff^r) = r, x^„(t) = 0(rs), we get ^gAsxj^ if^ = e + 0(r?), namely that to this 
order r is the proper time of any observer of the congruence. 

Yet, the ADM canonical formalism gives us an additional information. Actually, on each 
space-like hyper-surface Y,^^^^ of the foliation, there is a privileged contravariant space-like 
direction identified by the lapse and shift gauge variables 

Af^{T,a) = -^^-i— (0;n'-(r,a)) = (O; n.(r, a)) + O(r^), 
\N{T,a)\ ^ ^ ^ ^ 

AfA{r,a) = \N{r,a)\ (l; ^^^) = (|n(r, a)|; n.(r, a)) + O(r^), 
N\t, a) hir, a) = 0, Af^ir, a) N^{t, a) = -e. 



\N{T,a)\ = ^{^grsNrN^){r,a) = \n{T,a)\ + 0(4). (7.15) 

If 4-coordinates exist, corresponding to an on-shell complete Hamiltonian gauge fixing, 
such that the vector field identified by A/" (r, cr) on each E(^^^) is surface-forming (zero 
vorticity ^^), then each T,r can be foliated with 2-surfaces, and the 3+1 splitting of space- 
time becomes a (2-|-l)-|-l splitting corresponding to the 2-1-2 splittings studied by Stachel 
and d'Inverno [42]. 

We have therefore a natural candidate for one of the three spatial triads of each observer: 
El^^^^r) = Miir) = U\TjsSr))- By means of /4 (r) = /^(t,p^„(t)) and AA/„(r), we 
can construct two null vectors at each space-time point 



^^This requires that MAda^ is a closed 1-form, namely that we have drrir = dr\n\ + O(r^) and 
dr Us = dgUr + 0{r'^). In turn, this requires rir = dr f + 0{r'^) with dr f = \n\ + const. Our gauge 
has a non-surface forming J\f^. 
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^iir) = ^ {liir) -^f,^^r)) = -L (/4 (r) - (r)) + O(r^). (7.16) 

and then we can arrive at a null tetrad of the type used in the Newman-Penrose formahsm 

[43]. The last two axes of the spatial triad can be chosen as two space-like circular complex 
polarization vectors E^^(^_^^{t), like in electromagnetism. They are built starting from the 
transverse helicity polarization vectors -^a^^ (i 2) (''")> which are the first and second columns 

o 

of the standard Wigner helicity boost generating ICg^ (r) from the reference vector JCt^ (r) — 
\N\ (l;00l). 

Let us call J^f ^^(r) the ADM tetrad formed by /^^(r), AfUr), E^^ (1,2) (^)- This tetrad 
will not be in general Fermi- Walker transported along the world-line a;^^(r) of the observer. 

Another possible (but only on-shell) choice of the spatial triad together with the unit 
normal to is the local WSW (on-shell) compass of inertia defined in Eqs.(7.4). This 

local compass corresponds to the standard of non rotation with respect to the fixed stars and 
its T-evolution, dictated by Einstein's equations, does not correspond to the FW transport, 
which is defined independently from them using only local geometrical and group-theoretical 
concepts. The WSW local compass of inertia corresponds to pointing to the fixed stars with 
a telescope and is needed in a satellite like Gravity Probe B to detect the frame dragging 
(or gravito-magnetic Lense-Thirring effect) of the inertial frames by means of the rotation 
with respect to it of a FW transported gyroscope. 

Given the 4- velocity ^^(t) = Eg^{T) of the observer, the observer spatial triads Eg^ (r), 
a = 1, 2, 3, have to be chosen in a conventional way, namely by means of a conventional as- 
signment of an origin for the local measurements of rotations. Usually, the choice corresponds 
to Fermi- Walker (FW) transported {gyroscope-type transport, non-rotating observer) tetrads 
E^ffl^J)^ {t) 1 such that (we show also the implication of the linearization) 

C'-'^'M = 4.(t) ll(r) - al(T) li(r). (7.17) 

The triad Eg^^'^J^'^ij) is the correct relativistic generalization of global Galilean non-rotating 
frames. Naturally any other choice of the triads (Lie transport, co-rotating-FW transport,...) 
is possible. A generic triad ^#„(a)(r) wiU satisfy ;^^#„(a)(r) = ^^^/^(r) E£(„)(r) with 

^ ^fW)AB ^ ^iSR)AB ^.^^ ^^^^.^^ ^^^^^.^^ p^^^ ^iSR)AB ^ ^ABCD ^^^^ j^^ 

Jg^ Iffo A — 0, producing a rotation of the gyroscope in the local space-like 2-plane orthogonal 
to li and Jg . 

In the linearized theory on the WSW hyper-surfaces of our gauge FW transport implies 
no r-dependence for the triads -E^f (^)^'^(t). Neither the ADM tetrads nor the WSW tetrads 
(7.4) are FW transported. 

See Ref. [41] for the description of a geodesies ^^(t), the world- line of a scalar test 
particle, from the point of view of those observers ^ao,y{T) in the congruence which intersect 
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it, namely such that at r it holds x^t^ y{T)i^) ~ V^i'^)- The family of these observers is called 
a relative observer world 2-sheet in Ref. [41]. 



C. The Null Geodesies, the Deformed Light-Cone and the Eikonal. 

The solution of Eqs.(6.1) for null geodesies has still the form 
t{s) = + + s - dsi £ ds^'r^a + b^o)S2)b'b^ 



a-{s) ^a^ + b-- ds^ £ ds2 - 2 ^t^^^ia + 6(o)S2)6[o) b^^ - ^f;f,(a + b^o)S2Wb' 



with 



dT{s) 



ds 

da^('i) /■« r - - 1 

= = fe" - I ci52[(&[o))'-2'r«,(a + 6(„)S2)&[.)6'--'r;f,(a + 6(,)S2)^^^^^ (7.18) 



bU + bl) - [ ds2'tl^{a + b^b\ 



but now at s = we have 



^gABP'^P^\s=0 = e 



{p^{s)y - 2nr{s,ra{s))f{s)p^{s) - [Srt - e%rt{s,ra{s))f{s)p\s)] 



\s=0 



= 0, 



= (6^)2 - 2¥¥nr{a\ r-a{a^)) - [5rt - e%t{a\ rs(a^))] b''b\ ats^O, (7.19) 
whose solution is 

b^ = + == b'riria-', r,(aO) ± ^{b-nr{a-, r,(a-)))2 + [5,, - e^KM, r,(a-))] 6^6*, 

(7.20) 



= n.(a^ r,(a^)) 5^ ^ ? '"-^"^'M"^,^^^^^^ + 



2 

Therefore, if we consider the family of null geodesies emanating from a fixed point of 
space-time and parametrized by {b^}, we obtain the null surface describing the deformed 
light- cone through that point in our post-Minkowskian Einstein space-time. 

Let us now consider the eikonal equation [7] ^g^^{r, a) Da U (r, a) ds U (r, a) = 0, whose 
solution, the so-called optical function U, is used to find the null hyper-surfaces U (r, a) — 
tangent to the deformed light-cones, generahzing the planes — ± x^) — tangent 

to the light-cone in special relativity, where they are used for the front (or null) form of 
dynamics in light-cone coordinates. 
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In Minkowski space-time with Cartesian coordinates the eikonal equation is 
{drUM^T, a)y — {dUM{j,(y)Y = 0) namely drUu — Q;|5C/m| with a — ±. The solutions 
are Uuir, a) = /(r ± Er A a^) with = 1. 

After the linearization in our gauge the eikonal equation becomes (7''* = 7^ S^^) 
(dr Uir, a)f - 2 dr Uir, a) ^ n,(T, a) dr U{t, a) - ^ 7,(t, a) (9, C/(t, a))^ = 0. (7.21) 

r r 

Let us put U — Um + V with F = 0{ra). By disregarding terms of order V'^ the eikonal 
equation becomes the following quasi-linear partial differential equation for V 

2\^a\dUM\—'^nrdrUM drV — 2'^\a\dUM\'nr + '^rdrUm drV — 

r r 

= E (7. - 1) (dr Um? + 2a\dUM\Y.^rdr Um. (7.22) 

r r 

If we write this equation in the form a dr V^+Er drdrV = F, then a solution U = Um + V 
of the linearized eikonal equation can be obtained with the method of characteristics [25] if 
it is possible to find an exphcit solution of the system 

dV dr d(7^ ,„ 
— = — = . (7.23) 

F a Gr 
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VIII. CONCLUSIONS. 



We have defined a background-independent Hamiltonian linearization of vacuum canon- 
ical tetrad gravity in a completely fixed 3-ortliogonal gauge in the framework of its rest- 
frame instant form, where the evolution is governed by the weak ADM energy. In this 
non-harmonic gauge every quantity is expressed in terms of the DO's ra(r, a), 7ra(r, a), 
a = 1,2, parametrizing the independent degrees of freedom of the gravitational field. The 
method is based on the linearization of the Lichnerowicz equation for the conformal factor 
of the 3-metric and on the restriction of the weak ADM energy to its part quadratic in the 
DO's. 

As a consequence we succeed for the first time in solving all the constraints of tetrad grav- 
ity (super- hamiltonian constraint included) and to find a solution of the linearized Einstein 
equations in the uniquely defined 4-coordinate system induced by the chosen gauge, which 
corresponds to a post-Minkowskian Einstein space-time of the Christodoulou-Klainermann 
type. The DO's ra(r, a) turn out to satisfy the massless wave equation even if we are not 
in a harmonic gauge and we get non-zero shift functions (namely our 4-coordinates are 
non-synchronous). Besides re-opening the Hamiltonian approach to gravity, we can show 
explicitly the role played by the deterministically predictable DO's of the gravitational field 
in deforming the structures of the fiat Minkowski space-time without having used it as a 
background for the propagation of a massless spin 2 field as it happens in the standard 
treatment of gravitational waves in the harmonic gauges. In particular the two configu- 
ration DO's ra(r, a), a = 1,2, of our gauge, where the 3-metric is diagonal, replace the 
two polarizations of the TT harmonic gauge (with its non-diagonal 3-mctric) and, through 
the geodesic deviation equation, they induce two well defined patterns of deformation on a 
sphere of test particles. Even if we were able to eliminate the background with our Hamilto- 
nian linearization, our results are still coordinate (i.e. gauge)-dependent like in all existing 
treatments of gravitational waves: indeed our DO's r^, tt^ are coordinate-dependent non- 
tcnsorial quantities. To get a coordinate-independent description of them, we have to verify 
the main conjecture of Ref. [15]. It states that by means of a Hamiltonian re- formulation of 
the Newman-Penrose formalism [43] it should be possible to find a Shanmugadhsan canon- 
ical basis in which the DO's are also (coordinate-independent) Bergmann observables and 
also the gauge variables are coordinate-independent. 

Wc have also made some comments on the coordinate dependence of the gravito-magnetic 
effects: for instance in our non-harmonic gauge the gravito-electric-magnetic analogy does 
not hold. Again the verification of the main conjecture of Ref. [15] would allow a coordinate- 
independent description of gravito-magnetism. 

These results about background-independent gravitational waves in post-Minkowskian 
space-times are welcome because they open the possibility, after the introduction of matter, 
to study the emission of gravitational waves from relativistic sources without any kind of 
post- Newtonian approximation. For instance this is the case for the relativistic motion (but 
still in the weak field regime) of the binaries before the beginning of the final inspiral phase: 
it is known that in this phase the post-Newtonian approximation does not work and that , till 
now, only numerical gravity may help. In a future paper [44] we will add a relativistic perfect 
fiuid, described by suitable Lagrangian [45] or Eulerian [46] variables, to tetrad gravity, we 
will define a Hamiltonian linearization of the system in the completely fixed 3-orthogonal 
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gauge, we will find the Hamilton equations for the DO's both of the gravitational field and of 
the fluid, we will find the relativistic version of the Newton and gravito- magnetic action-at- 
a-distance potentials acting inside the fluid present in the weak ADM energy and finally, by 
using a multipolar expansion, we will find the relativistic counterpart of the post-Newtonian 
quadrupole emission formula. 

Let us remark that till now we have a treatment of the generation of gravitational 
waves from a compact localized source of size R and mean internal velocity v only [22] for 
nearly Newtonian slow motion sources for which v << c, >> R: outgoing gravitational 
waves are observed in the radiation zone (far field, r >> ^), while deep in the near zone 
{R < r « ^), for example r < 1000 ^, vacuum Newtonian gravitation theory is vahd. On 
the contrary with our approach in suitable 4-coordinates we are going to obtain a weak field 
approximation but with fast relativistic motion in the source subject to the restriction that the 
total invariant mass and the mass currents are compatible with the weak field requirement. 
This is enough to get relativistic results conceptually equivalent to the re-summation of the 
post-Newtonian expansion. 

Moreover we will have to explore whether our Hamiltonian approach is suitable for doing 
Hamiltonian numerical gravity on the full non-linearized theory. 
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APPENDIX A: RESULTS ON TETRAD GRAVITY. 

In this Appendix we reproduce those results of Ref. [3], which are needed in this paper. 
The Lichnerowicz equation in the 3-orthogonal gauge 7t^{t, a) ~ given in Eq.(193) of 
Ref. [3] is (see Eq.(A4) for the definition of the T^^y,'s) 



1.^ 



U 5 

J r I 

E '^ia)r (B, B„ t) E TS^TT-,) (r, a,) x 

b 

E ^2, r) (0-2e-75 ^ 7,,7r,) (r, a^) + 

s c 

b 



(Al) 



This equation for the conformal factor 0(r, a) of the 3-metric is imphed by the super- 
hamiltonian constraint, Eq.(191) of Ref. [3], which is 



T-Cnir, a) |^ = e0(i", B) 



67rG 



+ 



X 



/ dVi^5f,)7^«),(a,ai,r)(r'e-^^^^^'^'^''EM)(r,aO + 

r I 

+ r'{r,B) I dVidV2(Ee^^^™^"'"^ X 

^ b 
E ^2, r) (0-^6-71 ^ ^ 

s c 

I „~75 (7au+7a-u)f'a(T,CT) / r« X" X^" N 
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b 

E '?2, r) ^"-^ E 7c-.7rc-) (r, ^2) ) ] ] «^ 0. (A2) 

s c 

The scalar curvature and the Laplace-Beltrami operator appearing in these equations, 
Eq.(190) of Ref. [3], are 

+ 



^R[(f), Ta] = - E{(29«^^ <^ + ^ E laudvra) {4:djn - ^ E Tbudvn) 
2 1 



1 2 
+(29„Zn + ^ E lauduTa) C^djn (p ^ E Tbudun)] , 

r V<3 5 

If 7r(^)(T, ~ and 7f^(r, a) ~ are the Abelianized momenta corresponding to the 
rotation and super-momentum constraints, the connection of the old cotriad momenta to 
the final ones, see Eq.(1.4), after the quasi- Shanmugadhasan canonical transformation (1.3), 
i.e. Eq.(185) of Ref . [3], is 



= E/^'^i^(a).(^'^i;^)n^(^'^i) = 

s 

s 

[^ + V^E7-6.7r-J (r,ai), 

b 

with the kernel 
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(A4) 

The canonicity of the transformation (1.3) imphes that the kernels ^o)s '^ determined by 
the kernels G(a)'s. 

Eqs.(156) of Ref. [3], giving the partial differential equations for the kernels G^^j's, are 



g«(r,a) aq„")(ai,a;r) 



1 dG'^{d^,a-T) 1 9G^:5(ai,a;r) 



+ 



(a)' 



+ 



"(a)r2"ri 



(r, a) QrAr,a)Q 



1) s = a homogeneous equations : 

1 9Gg)(ai,a;r) 1 g; r) _ 

Ql{r,a) da^ Ql{T,a) da^ 

1 9Gg)(a„a;r) 1 aGg)(ai,a;r) _ 

Qi(T,a) Qi(r,a) 

1 dGll^(a„a;T) 1 aGg)(ai, a; r) _ ^ 

Qi(r,a) da^ Ql{r,a) da^ ° 

2) s 7^ a [s 7^ r, r 7^ a] in-homogeneous equations : 

1 aq:)(ai,a;T) 1 a; r) _ 

Q2(r,a) dar Ql{T,a) da- 

1 aq-)(ai,a;r) 1 aq:)(ai,a;r) _ ^3(^^.g) 

Eqs.(159) of Ref. [3] for the kernels /C(a)u's, which ensure the simultaneous satisfaction 
of the rotation and super-momentum constraints, are [Qu = 0^ ev3 7auraj 
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r 



+ 



+E 



= 0, a^h, 

= (<5(a)(6)f^r + e(a)(fe)(c) ^tUr(c) (t, o")) /C(5)„ (a, ; t) 

duQr{r, a) 



= ((^(a)(6)^r- + ~ '^(a)^'^(b)r-) !^ L U ) 

duQrir, B) 



-<5«(t,(7i)- 



0, 



a) (''■) ^) = X] ^(o)M(n)(^(m)r-(^(n) 



(A6) 



The cotriads, the triads, the 3-metric, the Christoffel symbols and the spin connection 
on the WSW hyper-surfaces of our gauge, given in Eq.(184) of Ref. [3], have the following 
expressions in the canonical basis (1.3) 







e(a) 




3 

9rs 













3f^r 



-5uv E '^le^ [2djn + ^ E Tbu^.r-, 



+ 



v3 5 

^ra->0 e(a)(6)(c) 5(6)r'^(c)« 25„/n </). 



(A7) 



Eqs.(102) and (105) of Ref. [3] define the Green function of the covariant divergence 
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m^)r{r. CtU^. B'- t) = B'). (A8) 



(a)(6) 



= 5(a)(6) c'' (a- ai), (A9) 
since the linearization reduces d!^{a,ai) to the Green function of the ordinary divergence 

The weak and strong ADM Poincare' generators of Eqs. (25) of Ref . [3] are 

^3 



ADM = J d ^4y^^v7 9 ( r^e r^n - r^.- l^n) 



S"""^ 3ri ^firs 3fvuv]f^ -i\ 
Grsui; ii ii \[T,(T), 



P'adm = -2/ d'a'VUr.BYW^T^a), 
J ADM — ~Jadm — J d^o-€{a^ 

Padm — Padm + j d^'7'H{T, a) ^ Padm^ 
Padm = P'adm + y d^a^'H^{T, a) ^ Padmi 
Ja^dm — Ja^dm 2 / d^'^'^^ 'HIt, a) Jadmj 

J ADM = J ADM + / d'a[a' 'n'{T, a) - a' 'H^t, a)] ^ J%^. (AlO) 

Eq.(197) of Ref. [3] identifies the Hamiltonian of the rest- frame instant form of tetrad 
gravity with the weak ADM energy 

-^(d)1dm = -^Padm,r = Eadm- (All) 

The expression of the weak ADM Poincare' generators after the quasi-Shanmugadhasan 
transformation (1.3) and the restriction to the 7r<^(T, a) = 3-orthogonal gauge, given in 
Eqs. (227) of Appendix B of Ref. [3], is 



Padm,r = e / ^V(y^['/>^E^~^^^" X 
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4 z 
- -j^drln(t)Y,rhrdrri + -{Y^TbAn)' 

^ b b 

^ a 



X 



+ 



" b 

I d'a, E 5ra)^a)r(^' ^1' ^l'^' ^^^^^ E TS^vr-,) (r, a 



^l) + 



6 

+ 



s c 

E %.(^, ^1, r|0, rj (0-^6-^ ^-^'^''^ Y^Tbrn) (r, ^i) 
^ b 

E ^2, r|0, rj (0-^e-7f S.^-- ^ ^..tt,) (r, a,))]), 

s c 

b 

V3 5 

E ^crTTc] (r, ct) + 

c 

+ V3 / dVi E [ - E (e^^^^^''"-^^-)'-n29,Zn0 + -^T.rbudrn))(r., 

•' s u V <-> h 



+ 



UV V <J f. 



+ 5;(29Jn + ^ E Tbvdun))){r, a) 



(0-V^^^^--E7c-.7re-)(r,aO), 
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J' 



b 

u V>J 5 c 



c 

b 

- A^9r(l)Y,rbrdrn+^{Y,rbrdrr-bf)]{T,a) - 

-^'^"V,^){2r^E^D(^.^) + 



+ 



+4 



a) X 



+ 



s c 

+ E^^ i\bf(a) - 0(a)0{b)) X 

uv 



X 



+ 
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m J, 



ec 

8^ 

2_ 

e V3 



^/3 Sa^''''''' To")''- ' ' "'■ 



{djn(f)+ ^JZiTbv- rbu)dur-b)\ (t, a). 



(A12) 



Eqs.(194) of Ref. [3] give the equation for the determination of the lapse function n{T, a) 
i consequence of the time constancy of the gauge fixing p(r, a) = ~ 

drp{T, a) = {p(r, a), H(d)adm,r} = J (fainij, ai){p(r, a),nR{T, ai)} + 
+ K{r){p{T, a), PXdm,r} + K{r){p{r, a), Padm,r} ~ 



+ 



_l_ ^ ^Padm,r I ^ (^-^ ^Padm,r 



+ \r{T)- 



S(l){T,a) ' ' 50(T,a) 
n{T, a) - n{T, a\ra, tt^, Xa] ~ 0, 



0, 



dr [n{T, a) - n{T, a\ra, tt^, A^]^ 

= A„(t, a) - {n{r, a\ra, tt^, Xa], H(j^fJ^M,R} ~ 0, 

=^ A„(t, a) determined; 



the rest-frame instant form expression of this equation is 



r, (Ti) ^Padm,r 

y^"^ "^"'"^^ c50(r,a) 



(A13) 



For 7rtf>{T,a) — the shift functions nr{T,a) of our gauge, as functions of the lapse 
function n(T, ct), have the form given in Eq.(187) of Ref. [3], i.e. 



;4V37rG 



(r, a) / o?^ai[e -n(r, ai)]0 ^(r, ai) 



tow 



(A14) 
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For TT^ = the extrinsic curvature of our WSW hypcr-surfaces, the ADM metric mo- 
mentum and the DeWitt supermetric, given by Eqs.(186) of Ref. [3], are 

3i^,,(r, a) = ii^[e75 E.(7c-.+7..)r-. ^(^^^^^^^^ + ^^j^^^^ _ ^^^^^^^ j 

^ u 



2 1 ^ o 

3fV _ r3-r S-'* 



%/3 



(a) 



E / (iVi/C[„)„(CT,<Ti,T|(/),rs,n] + 



+ e 



6 

/'_;5N_r3;^ 3- _i_3- 3- _3- 3- \(^p±\_ 
^rsuvy'i^j L yi'u ijsv ' ijrv ifsu ifrs yuv\\'i'^) 



(A15) 



The tetrads and cotetrads adapted to the WSW hyper-surfaces of our gauge, given in 
Eq.(l) of Ref. [3], are 



4 p^t 
(S)^(a) 


= iUK) = ^' = 


= ^(&^ 




4 p(a) 




= K = Nbl- 




4 E^M 

(S)^(a) 








(S)^(a) 


-(S)^(a)^' ^ 


4 

(S)^(o) = 


A 
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(S)^(o) 



(S)^(a) 



-(»)' 



4 ~(") 



4 



x{o) 



4 



4 



x(a) 



3p(a) 



(A16) 



The 3-metric on our WSW hyper-surfaces and the hne element of our space-time, given 
by Eqs.(183) of Ref . [3], are 











— T 



0' '^'r^ 



1 ^ 
q = 2ln(i) = - In g, 
6 



'9 



(A17) 



The Frauendiener equations for the WSW triads, given in Eq.(200) of Ref. [3], are 



3 {WSW)r 



(1) 



3 (WSW)r 



3 {WSW)r 



r ^(2) 
3i 



= 0, 



(3) 



3 (WSW)r 3 (WSW)s Svy 3AWSW) 3jWSW)r 3(WSW)s Svy 3JWSW) 
^(1) ^^(3) 6(2)3 '^(S) ^(2) ^ '^(l)s 



3^(W5W)r 3^(l^5U^)3 3V7 3^(W^SVF) _ ^ 
-h 6(2) e^^) Vr 6(3)^ — U. 

The WSW tetrads of Eq.(201) of Ref. [3] are 

1 



(A18) 



4 



-e -|- n 



4 



(a) - lU, e („)j. 
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UK)-b's (A19) 

The embedding of the WSW hyper-surfaces of oue gauge into the space-time, given in 
Eq.(208) of Ref. [3], are 

^^wsw)ir: B) = 5f,)x[£)(0) + P(r, B)t + ejf (r, B)a^ = 
= ^foo)(0)+&!^(r,a)F^(r,a), 



-e + n(r, a) 



<f ) = 'eC-^-'f., (T, - Z^^J^T, (A20) 

with the transition coefficients b'X determined by Eqs.(209) of Ref. [3] 

or y^=(A-')B''F^^. (A21) 
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APPENDIX B: THE KERNELS F^^^^y 

As shown in Ref. [3] the kernels Fj^allb)^^ curvature hmit of the kernels F^'^^^^^'s of 
Eq.(1.4), are given by the equations 



Ha){s){b)0 (CT, CTlJ - e(«){r)(6) 

u,r ^^1 



-,(o)r X ^^(a) (^1'^;''") . rrrS/- -\ 
'(a)(l)('^l'^;^) = ha)^^2^ i^^^l) = 



= ()(„)i()2<) (cx,ai) ^-^^2 



r 



so that from Eqs.(2.12) we get 



a, 6 = 1, 2, 3, 



'(3X2) (^'-^l) = -^(1X2) (^^-^l) = 



= ^(1X3) (^'-^l) = -F{2)(3)i^^^l) = 



^(ixV^'-^l) = ^(2Xl)('?''?l) = 

1 Q 



1 (9 /■'^i 
= "2^K'^ )^(^ '^i) ' 



^(2X3) (^'^l) = ^(3X3) (^'^l) = 
1 a 



-2 5^71'"^''^"^'"^'"^"^) = 



^(ixV^'^i) = ^(ixV^'^i) = 



(Bl) 
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^(o)3 
(2)(2) 



^(o)2 
(1)(2) 



^(lK3)(^''?l) 



1 d 



2 dai J-. 



1^2 3\ 



2^^^!'*^ ) 0{O- ,(7j, 



^(0)2 
(3) (2) 

1 d 



1 d 

2 
1 



2 2^'95(a^cT?) 



1 a 



2 9(7? 



"1 
00 



1^2„,,3\ 



1 1 iNl^fK^^-^l)/!/ 3 3^ 
-^^(^S^l) g^2 ^ (^^^ )• 



(B2) 
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APPENDIX C: FOURIER TRANSFORMS. 



By assuming the validity of the Fourier transform on the hnearized WSW CMC-hyper- 
surfaces, which are assimilated to flat surfaces at the lowest level, for the real functions 
ra{r, a), TTair, a) we have 

r-a{r, a)^ j -^^f-a{T, k)e'^-^, f?(T, k) = fa{r, -k), 

7ra(r, a) = j ^^7fs(r, k)e'''-'^, 7r*(T, k) - n-a{T, -k), 

f,(T, k)^l d'ar-a{T, a)e-''-'', 7r,(T, k) ^ J ciV7r,(T, a)e-'^-^, 

= 5^1 J dVi J dVe-*-^ie-^^-^5(a - <ti) 

= 1 d'ae-^^^+^'>' = {27r)%iS{k + k,). (CI) 

The conditions (2.15) and (2.16) become 

7ra(T, O/c^/c^) = 7ra(T, /c^O/c^) = 7fs(r, k^k'^Q) = 0, 



fair, Oee) = fair, k'oe) = fair, k'eO) = 0. (C2) 



Some useful relations are 

f dk e*''^ f d^k e'^'" _ 1 

V'(x) = y dy0{x - y)f{y) f{k) = iki^{k), 

^''^^^'^^ = / f^t^' + |^|>~.(r,^)e^^^-. 



/(c?i) /(^) 



i J 47r|(7-(7i| |A;|2 

J Unla-aA ~ J (27r) 



47r|a-ai| J (27r)3 |^|2 
/ da^,eia\ a^) / ^^^^(a^", ar)e^^"^^" = / da[^(<7[, <j")e^^^'^^ = 



Since 
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' ' ' doll dale 



oo J — oo 



CO fOO 



dal / da^ / dal / da^ / rfVs 



^lr^i^L^3Vc-(-r, ^3) 



47r|(T2 - <?3| 



/ 



(C3) 



Let us consider Eqs.(5.2). Their Fourier transform is 

drfa{T, k) = ^A^-^{k)7ri{T, k), 



(C4) 



where the matrix ^^^(A;) is 



Aab{k) = Al^{k) = A^l{-k) = — S^l+- ^ 7ar76. 



r,s,u,v 



det (A^i 



c6 {k^k^k^y 



(C5) 



Eq.(C2) imphes that Eq.(C4) is well defined, even if the matrix A^iif^) diverges for 
k^ — > 0, if 7ra(T, k) vanishes for k ^ at least as {k^k'^k^y'^^, e > 0. 

For A; 7^ its inverse is 



detA{k) 



cd 



1 

SttG 1^1 6 



{k'k'e)%l + 



+ 2 E e-a^-^eulds{'i5tr-l){25ts-lWk'ky 

r,s,t,c,d 



(C6) 



The Fourier transform of Eq.(5.5) is 

drTfa{r,k) =Y,Bab{k)fi{T,k), 



(C7) 



— * 

where Bg^{k) is the matrix 
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satisfying 



det (b^i 



T rs I "'I 



{k^k^k^'f 
1927r2G'2|^p 



Y.A,j,{k)B-,,{k) = -\k\H-a, = J2B,i{k)A-,,{k), 



(C8) 



A-\k) = -\k\^B{k). 



(C9) 



But this implies that the Fourier transform of Eq.(5.7) reduces to the wave equation 



jj r rs \K\ 



c rs tw uv 



[l-2{5ur + 5,rm-2{6^s + S,s)] 



h.2h.2 l2l2 



3' 

)^c(t, A;) 



= -|^|V~,(t,^), 
whose solutions 



(CIO) 



d^k 



(27r)3 



(Cll) 



have the arbitrary functions Ca{k) satisfying Eq.(C2). Then Eqs.(5.4) imply the following 
form of the momenta restricted to the solutions 



TTair, a) = -i 



d^k 
(2^ 



Ur,k)^j:A;iik)drhir,k)^ 

b 

= E^a^ (^) \ - + t\k\C-;{-k)e^\'\^ 



(C12) 



To have TTa{r, k) vanishing for A; ^ at least as (k^k'^k^)'^'^'^, e > 0, we must require that the 
functions Ca{k) vanish at least as {k^ k"^ k^Y for k^ — > 0. 



The Fourier transform of the ADM energy (3.1) is 
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Eadm = / T^-T^ H^Tslr, k)7rair, -k) + 



l^'^'J ^j; r,s,u,v 



[1 -2(5„, + — J 7fa(r, A;)7rb(r, -/c) - 



247rG / T^X'^--^~brk''Mr, k)n{r, -k) + 

^ ^ abr 

+ / ^J--r^^'-'^r, k)nir, -k). (C13) 



As a check we can recover the Fourier transforms (C4) and (C7) of the Hamilton equations 
by using this form of the weak ADM energy and the Poisson brackets (CI). 

The Fourier transform of the shift functions (3.5) is 

OttG k k^ 

nr{r, k) = _E 75.(1 - 5^,)[l - 2(5„, + 4.)][1 - 2((5„, + 4,)]^^5(r, k), (C14) 

hvuc ^ 

while the Fourier transform of the extrinsic curvature is 

^~krs{T, k) = ev^^ [^^rsTas - (1 - 5rs)Y.Taw[l - + (^.«,)] 7^] ^a(r, k), (C15) 

The following Fourier transforms are used in the study of the geodesic deviation equation 

'X ^ — * A — * 

^r^^(T, k) = iksUrir, k) + e^Krs{T, k), 



%,At, k) = -edr'Krsir, k). (C16) 



The Fourier transforms of the exponent q{T,cx) (v^ ~ ~^(1 + Sf/)) of the conformal 
factor of the 3-metric, given in Eq.(2.2), and of the inverse 3-metric (3.7) are 

q{r, k) = ^ E7.n^r,(r, k), (C17) 

and 

T{r,k) = -^Y.b-arfa{r,k) + lY.TauS-/,{T,k)]) . (C18) 
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APPENDIX D: LINEARIZED 3- AND 4-TENSORS. 



1. Linearized 3-Tensors on the WSW Hyper-Surfaces E,- 



/,Prom Eqs.(3.8) and (A7) we get the following results for the hnearized Riemannian 
structure of the WSW hyper-surfaces 
The Christoffel symbols are 



^fL(r,a) = i=X:(-5, 



1 



laudrT-air, a) - 



J 47r|c7 — (7i| 

J 47r|(7-(7i| 

+ ^rv [lavduTair, B) - 

J 47r|(7 — (7i| 



+ 



+ 



)+0(r^), 



V ^ aw -I 47r|(7-(7i| 



47r|CT — (Ti 

E ^r^Jr, ^) = 4^ E ([ - ^uv(l-au^ + ^ ^ 7a«;5y + (7sn + Tav)dudv 



»^a(r, a) - 



i 47r (J 



47r|a — ai| 



-ilV Vtl ^ 111 



(Dl) 



while the spin connection and the field strength are 



ia){r,B) = ^ Ee(a)(r)(n) E [larduTair, B) - 



47r|(T — (Ti| 



(D2) 



^^^rs(a)(T, o") = (9r ^'i^sia){r, B) - 8^ ^a'r(a)(T, B) + 0(r?) 

= 4? E [lasdudrTair, B) - 
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J 47r (7 - (7i 



'(a)(r)(u) 



J 47r|(7-(7i| 



0(r?)- 



(D3) 



Finally the Riemann and Ricci tensors and the curvature scalars are 

^Rrsuvi''' , B) = e(r)(s)(a) ^^ui)(o) , B) + 0(r?) = 

= X] ([7ar((^rt;9s9„ - Srudgdy) - ^asi^sydrdu - Sgudrdy) ra{T, a) 

2 ^ ^ Tat [^^rt;^s^M ^ru^sdy (Sgydrdu ^su^^r^^v) 
via ^ t 



r V "J 



(T,a) + 0(r^) 



(D4) 



2. The 4-Christofrel Symbols, the 4-Riemann Tensor and Einstein Equations. 

In this Subsection we will give the linearized form of the main 4-tensors of our space-time. 
By using the paramctrization of the 4- metric ^Qab given in footnote 6 and Eqs.(4.1), 
with the ^r^s'^ of Eq.(Dl), we get the following 4-Christoffel symbols [A^ = — e + n, — 

rir + 0{ri)] 

^f;,(r, a) = {^[d^N + N^drN - N^N^ 'K^s]) (r, a) = 
'KAt: B) = {^[drN - 'krsN^]) (t, a) = + O(r^), 
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+ [1 - (^rs] X^7o«)[l - 2{5ryj + Ssw)] 

w 



3 f,uv 



N 



)Nd^N + 



+ N\N^ - 2N{T - ^^fK^rN^) (r, a) = 

^drnu{r,a) + 0{rl) = 
-27iG 



d rT" 



at ran 

da 

AT" 



/a fa roo roo 

da'C / da"^ / dcr^ / dcr^ 



= a,n,(r, a) + e^Kur{r, a) + 0(r2) = 



drTiuir , a) H ^^Go(a)(b)(c)(d)0( 



E [25"7.r7r.(r, a) + [1 - 6^] E 7.41 - 2(5," + ^rt)] 

^ a t 
fp' roo POO 



S(a)S{b)r E ^ic)w^^7d)i^^ B) + 0{rl) 



+ 



1 



+ ^ E 7.t E(l - <^-n) [1 - 2(C + [1 - mm + <5tn)] 

^ at mn 



da'da'- 
+ 0{rl), 



da"^ / ci< / da: 



da. 



{daif 



+ 



AT" 



df^rair, ai) 



47r|a — (Ti| 



+ 



+ ^ruhauOvra[r, CT) - T,l^lawOv / d (7i — — — 

L I'rf J 47r|(j-ai| 



+ 



+ 6Jl-a.dur-a{T,a) -W^-,^du I d'a, + 0{rl). 



47r|(7 — (7i 

Since (t,ct) = in our gauge, the space-time spin connection ^a;^*^"^/?) 



(D5) 
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^E^ldA ^£^5) + ^Ef^^] can be evaluated with the formula 



(D6) 



by using Eqs.(3.9) and (D5). 

For the 4-Riemann tensor ^R'^bcd = | ^9"^^ (^b^d ^Qec + dEdc ^Qbd - dsdo ^Qbc 
dsdc 'qed) + 'g^"^ '9EF {'Tie 'Tin - 'TId 'T'^c) we get 



4 5A 



*R\cD{r,a) = + O{rl), 
'R\rs{r, a) = dr 'tl, + 0{rl) = ed, 'Krs + 0{rl] 

= V^^- 5Z [2(^rs7ar7ra(r, ct) + 
+ [1 - 5„] ^7a^[l - 2(5^1^ + 6sw)] 



POO roo -| 



^ a t 
^ fi'^ roo roo 

^3 pco rco 



(r, w^'w^a'^^'''")) 



+ 0{rl), 

'R\sr{r, B) = -ed, 'K^sir, a) + O(r^) = 

+ [1 - (^rs] Xl7a«;[l - 2{6rw + ^sw)] 



a 



w\2 



/ dal daldr7r-a{T,a[aia''^^'^)]+0{rl), 

J a'" J 



9' f,3^ idlr,)ir,B) 



1 r „ „ , ^, 1^ /■ 3 



- S^TasdudrT-air^a) ' ^^Ta^g^ J ^ 47r|a-a,| 



+ 



4:7r\a — (Ti| 
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— Sr 



1 (9^ 



2^ "'"'da'' da' 



3^ (^lVa)(^,^) 



d'^ai 



+ 0{rl). (D7) 



47r|(T — a\\ 

Then we get the following linearization of the Ricci tensor and of the curvature scalar 



AB 



'R 



4 T-tE 



E ^it:V,. = -eE^r 'i^.. + O(rD=0 + 0{rl), 

r r 
r r 

eY,dr'Kru + Oirl)^0 + O{rl), 

r 

4 At _|_ \ ^ 4 fyu 

rrs ~r / ^ -ft ^ws — 

u 

u 

R = ^77^^ 'Rab + 0(rD = e(^^.. - ^ ^4.) + 0{tI) = 

r 

= -^E(^"'r;f, - a.3f;fj + o(r^) = o + o{rl). 



(D8) 



We see that the Einstein's equations ^.R^a=0 (corresponding to the super-hamiltonian 
and super-momentum constraints) are satisfied as a consequence of the results of Section II. 



The spatial Einstein equations ^Rrs = are independent from the shift functions rir and, 
after having used the Hamilton equations (5.6) to eliminate dr'Kaij^a), it can be checked 
with a long but straightforward calculation that they are satisfied 

'Rrs{T, a)^edr 'Krs + E(^« '^rs " d^s 'K)u){r, ^) + 0{rl) = 

u 

= V3^^ E \^5rslasdTT^a{r, a) + 
^ a 

+ (1 - 5rs)^law[^ - 2(5r«; + ^sw)] 

W 

dl da{ da{drT:-a{T,a{a{a''^''^') 



+ 



+ 4^ E ( f ~ ^rsilar^ + ^ E 7a«;<9^) + {lav + 7a5 )^r^J »^a(T, ff) 

- E lawdrds / d^ai 



dL^ajr, ai] 
Airla — (Ji\ 



^ aw J An a-ai 



(D9) 



Then we can obtain the field strength in our gauge by using the equation ^i^AS(a)(/3) 



4i7'4z7i"4d i34° i34° i4° 4°(')') 4° 40(1') 

(S)-^(a) (S)-^(/3) RcDAB^Oa B{a){l3) - Ob <^A(a)(/3) + <^A(a)(7) "^SC/?) " <^S(a)(7) '^A(^)- 
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3. The 4-Weyl and 4-Bel-Robinson Tensors and the Kretschmann Invariant. 



In this Section we will study the radiation gauge form of others 4-tcnsors. 
Let us evaluate the 4-Weyl tensor *Cabcd = *Rabcd + \{^Rac'^9bd - *Rbc'^9ad 
Rbd ^9 AC - ^Rad ^9bc) + 1{^9AC ^9bd - ^9ad ^9bc) ^R = ^Rabcd- We get 

^rstu -^rstu ~l~ ^( ^ -^rt^su ~l~ ^ -^st^ru ^ -^su^rt ~l~ ^ -^ru^st) 



— e- 



1 92 



,3, i9f^ra){r,a)i 



Ta 
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f rfVi 


(^L^a)(r, a) 


9(7*9(7'^ J 


4it\(j — (Ji\ 
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f d^ai 


idLr-a){r, a) 




47r (J — (Ji 


92 


1 dVi 


idLra){r, a) 




47r (T — (7i 



47r|CT — ai\ 
+ 



latdudsTaiT, ? ) " - $1 7a 

+ ^( - 1] [2(^rt7at^r7ra(T, ct) + 

+ (1 - X!7a«,[l - 2{5rw + 



+ 



V J 47r|(7 - (7i| ^ 



47r|(T — (Ti| 



+ 



^ a 

+ (1 - 5^t) ^7a«;[l - 2(5^^ + StJ] X 

/•oo roo -| 

Jcr^ Ja* 

+ 4^ II ([ - <^5t(7asA + ^ Il7a«;92) ^ ^ 7as) 9,9(1 ra(T, ct) - 



+ 



+ 



+ ^awdsdt j d^CTl 



47r|(T — (?! 



47r|(T — CTi| 
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+ 



+ (1 - 5su) X] 7o«)[l - 2((^s«) + Kw)] 

w 

+ 4^ H ( f - ^suilas^ + Tawdl) + in-as + lau)dsdu 

111 



47r|(T — (Ji| 



Vs. 



+ -?r y^lawOsOu 1 a CTi ^ + 

2 tlT 47r|a-ai| 
+ e(5stV^— ^ X] [2(^rs7a«9^7rs(T, ct) + 

^ a 

+ (1 - Sru)Ylaw[l - ^{Srw + ^uw)] 

w 

POO roo 

1 1 
V3 a ^ «; 

V 47r|(7 - (7i| ^ 

+ ^Z^Tawdrdu / — +C'(rJ, 

2 a«; 47r|(T-(7i| ^ 



ra{T, a) - 



-AttG 



a t 

/ 



1 

- (1 - 5ru)[l - 2{5rt + ^ut)]-g^^,y 

+ 0{rl), 



dw' (i«;"7rs(T, V^'"'" 



+ 



4^ _ 4p I _^4p 

^rrrs ^rrrs 2 



+ [1 - ^rs] X] 7at«[l - '^{Srw + ^sw)] 
w 

Qp' roo poo 



+ 
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+ 



+ (1 - Srs)^^aw[i - "^(Srw + 5sw)] 

w 

VO a ^ w 



ra{r, a) 



^ ^ J 47r cr — (Ji 



47r|a — ai\ 
3, df^M^^^i) 



47r a — ai\ 



0, 



(DIO) 



On the solutions of the Hamilton equations the Weyl tensor "^C^ BCDiT^ a) coincides with 
the Riemann tensor. 

Let us consider the congruence of time-like observers whose 4-velocities are given by the 

unit normals to T.^wsw)^ ^j^. \^^eI^-^{t,B) = el^{T,a) = -(1; a)), Js)^a V> ^) = 

Ia{t,^) = ~e(l;0) [sec Eqs.(2.2)]. The electric and magnetic parts of the Weyl tensor with 
respect to these time-like observers are 

^EAsir, a) = ^EBA{r, a) = ["Cacbd f 1'']{t, a), 



[^EaB l"] (r, a) = - \^EAr + ^EArUr] (t, g) = 0, '^^a(t, g) = 0, 

■^rsi.'^i [ ^TTST ~l~ ( ^rrsu ~l~ ^STru)'^u ~l~ ^rusv'^u'^vli.'^i 

— ^CrrsT + 0(r?), 

^Errir^B) = [^CrrTsnrns]{T, B) = - [^E^^^n^] (t, ct) = 
= + O(r2), 



^Hat + ^HatT^t 



^^*^C'yiBCD = ^CAB'^'^ ^C'efcd = f ecD^'^ ^C'abef is, due to the Lanczos identity, the unique 
Hodge dual of the Weyl tensor [47] (the Riemann tensor has different left and right duals) . 
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= + O(r^), 

= + O(r^). (Dll) 

Both '^Eab and ^^as have five independent components. See Ref. [47] for an analogous 
decomposition of the Riemann tensor. 

The four eigenvalues Aq, a = 1, .., 4, of the Weyl tensor are 

^ CD ^ EF ^ AB — ^[ra), 

'C^^^'^ ecDEF 'C'^^uv "C^'^AB = 0{rl). (D12) 

As shown in Ref. [15], the 4-coordinate system a"^ — {r, ct}, corresponding to our com- 
pletely fixed 3-orthogonal gauge, is identifiable by means of 4 gauge-fixing constraints 
a"^ — F^[Aa(r, ct)] ^ 0, where the are 4 suitable scalar functions of the 4 eigenval- 
ues of the Weyl tensor such that in the linearized theory we have F"^ = F^^ + 0(r?) with 
^(i)=0(r,). 

The Kretschmann invariant and pseudo-invariant are respectively [*'^R'^^'-^^ = 

I^ABEFAtd CDl 

-^e Kef \ 

A-A = 'RABCo'R^^'"'^0{rl), 

*R . R = ""R^scD = 0{rl). (D13) 

In Ref. [16] it is suggested that in presence of matter a coordinate-independent characteri- 
zation of gravito-magnetism is *R ■'B.^Q. 

The Bel-Robinson and Bel tensors are respectively 

in- — 4/^ E F I 4A E F 

■l-ABCD — ^AECF ^B D + '^AEDF ^B C ~ 

4A 4AEF G 4A 4A EFG , 

— 2 9AB '^EFCG <^ ^ ~ 2 '^AEFG + 

, ^4^ 4^ 4A 4(^EFGH _ n(r-'^\ 

o 

4}2 4£) 45 E F I 4 £) 4 £> E F 

£> — Kaecf j^b d + j^aedf j^b c — 

4p 4oEF G 4p 4p EFG , 

— 2 9^>l-B ^EFCG -n- ~ 2 ^^EFG ^B + 

J-^^n 4p 4t3EFGH _ (^l2\ 

+ g fl'AS S'CZ) i^EFGH ti —U\rg), 

^^ABCD = ^Tabcd + ''A^abcd + "^Qascd, (D14) 

where "^Qabcd = I'^R [^Cacbd+'^Cadbc] is the matter-gravity coupling tensor and "^TWascd 
is the pure matter gravitational super-energy tensor (see Ref. [47] for its expression). 
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